Broad relaxation spectrum and the field theory of glassy dynamics for pinned elastic 

systems 
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We study thermally activated, low temperature equilibrium dynamics of elastic systems pinned 
by disorder using one loop functional renormalization group (FRG). Through a series of increas- 
ingly complete approximations, we investigate how the field theory reveals the glassy nature of the 
dynamics, in particular divergent barriers and barrier distributions controling the spectrum of re- 
laxation times. First, we naively assume a single relaxation time for each wavevector k, leading 
to analytical expressions for equilibrium dynamical response and correlations. These exhibit two 
distinct scaling regimes (scaling variables Tk^ \nt and t/rk, respectively, with T the temperature, 6 

the energy fluctuation exponent and ~ e"'' ^'^) and are easily extended to quasi-equilibrium and 
aging regimes. A careful study of the dynamical operators encoding for fluctuations of the relaxation 
times shows that this flrst approach is unsatisfactory. A second stage of approximation including 
these fluctuations, based on a truncation of the dynamical effective action to a random friction 
model, yields a size (L) dependent log-normal distribution of relaxation times (effective barriers 
centered around and of fluctuations ^ L^^"^) and some procedure to estimate dynamical scaling 
functions. Finally, we study the full structure of the running dynamical effective action within the 
field theory. We find that relaxation time distributions are non-trivial (broad but not log-normal) 
and encoded in a closed hierarchy of FRG equations divided into levels p = 0, 1, . . . corresponding 
to vertices proportional to the p-th power of frequency u}^ . We show how each level p can be solved 
independently of higher ones, the lowest one (p = 0) comprising the statics. A thermal boundary 
layer ansatz (TBLA) appears as a consistent solution. It extends the one discovered in the statics 
which was shown to embody droplet thermal fiuctuations. Although perturbative control remains a 
challenge, the structure of the dynamical TBLA which encodes barrier distributions opens the way 
for deeper understanding of the field theory approach to glasses. 



I. INTRODUCTION 

Extremely slow dynamics is a ubiquitous property 
of complex and disordered materials. Despite many 
decades of research, current understanding of such glassy 
motion is Hmited to phenomenological models,^ mean- 
field theory and abstract caricatures in terms of 
the dynamics of small numbers of degrees of freedom 
in a complex energy landscapei^ Besides exact solu- 
tions in one dimension (e.g. for the random field Ising 
model^) little is known about the non-equilibrium be- 
haviour of realistic models. True disordered materials, 
from spin glasses to supercooled liquids to the pinned 
elastic medium, involve extensive numbers of local de- 
grees of freedom such as atoms and spins moving col- 
lectively in a random environment with either ex- 
ternal or self-induced randomness. The pinned elastic 
medium being the simplest model involving such physics 
we study it as a prototype. It is of interest by itself 
for numerous experimental systems such as vortex lat- 
tices in superconductors^S^ interfaces in magnetsjiSii^ 
charge density wavesji^ Wigner crystalsii^ The equation 
of motion is 



where u{r) is a height (or displacement) field, rj a bare 
frictional damping coefficient, c is the elastic modulus, 
f{u, r) is the quenched random pinning force, and (^(r, t) 
is a thermal noise. Here r is he d-dimensional internal co- 
ordinate of the elastic object. Both / and ( are Gaussian 
random variables with zero mean and second moment 



f{u, r)f{u', r') = A{u - u')5'\r - r') 
{ar.t)ar\t'))=2T^T5\r-r')6{t-t'), 



(2) 
(3) 



where T is the temperature and we set Boltzmann's con- 
stant kB — ^- The value of 77 generally sets the relaxation 
time scale, e.g. here rj — tocA^, Iq being the microscopic 
time scale and A the short scale momentum cutoff. In 
general, one may be interested in a variety of thermal 
and sample-to-sample fluctuations of the system, as well 
as various responses of the system to external probes. 
We will focus on the simplest of the latter, described by 
(linear) response functions. 



dUr{t) 



(4) 



in a given disorder realisation, and its disorder average. 



7]dtUrt = cVjurt + f{urt,r) + C,(r,t), 



(1) 



(5) 
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At equilibrium, both the single sample and the averaged 
response functions become time translationally invariant, 

"^1 — r'.tt' = T^i — r',t-t' and Rr~r',tt' = Rr-r' ,t-t' ■ 

The equatiuon |^ has the usual Langevin form, and 
guarantees the existence of a stable equilibrium prob- 
ability distribution, provided, as assumed here, that 
f{u,r) = —dV{u,r)/du is of gradient form. The equilib- 
rium distribution (strictly speaking defined in a finite size 
sample) has the Boltzmann form p{u) cx exp{~H[u\/T), 
with 



H[u] 



\Vu\'' + V{u{r),r) 



(6) 



Three universality classes of special interest are usually 
considered: (i) short range disorder A{u) which describes 
e.g. random bond (RB) disorder for magnetic domain 
walls (ii) long range disorder which describes e.g. ran- 
dom field (RF) disorder, and (ii) random periodic (RP) 
A(u) which describes pinned density waves or lattices. 
Although these systems differ in their details, e.g. in 
their roughness exponent u ~ r*" , they do not yield qual- 
itatively different behaviour in their dynamical response 
studied here. 

The aim of the present paper is to develop an ap- 
proach based on the renormalization group (RG) to study 
the low temperature dynamics of pinned elastic sys- 
tems described by Although we focus on equilib- 
rium dynamics, some of our considerations arc also rel- 
evant for non-equilibrium relaxation. It was shown that 
to describe the statics at equilibrium one needs to fol- 
low the full correlator of the random potential (or the 
random force) using a functional RG (FRG) method 
in a d = 4 — e expansion^S*ii Several extensions de- 
scribe correlated disorder,— the driven dynamics near 
depinningi2i2Si2i and the small applied force, thermally 
activated, creep regiment However until now the FRG 
has not been used to study the dynamical response and 
correlations in equilibrium or aging regimes, nor to probe 
the crucial question of the distribution of the relaxation 
times. These are important quantities directly probed 
in experiments where the system is often dominated by 
fluctuations or not able to reach equilibrium within the 
experimental time scales. We investigate this problem in 
three stages, of increasing accuracy (and, unfortunately, 
complexity), only in the last stage attempt is made to be 
exhaustive. A companion paper (Ref. Ila) is devoted to 
the statics. A short account of both works can be found 
in Ref. 24. Some of the present considerations concerning 
approximate schemes have also been discussed indepen- 
dently in Ref.^'"^ . 

The first question investigated in the present paper is 
the validity of the single time scale approximation within 
the RG. Specifically, in the first part of our study (Sec- 
tion II) we use, as was done in previous works fi2i22iS^ 
the simplifying hypothesis that the relaxation of each in- 
ternal mode, of wavevector fc, is controlled by a single 
relaxation time scale Tk- This allows to obtain closed 
equations, within the one loop FRG, for the general two 



time response and correlations, as measured in aging ex- 
periments. It yields, in the equilibrium regime on which 
focus from then on, interesting analytical expressions for 
the equilibrium response and correlation which exhibit 
two distinct scaling regimes with scaling variables In t 
and t/Tk, respectively (6* = li — 2 + 2^ is the energy fluctu- 
ation exponent, and ~ e'^*^ /^). However, several fea- 
tures of these results are found to be unsatisfactory, such 
as the non-monotonicity of the response as a function of 
wavevector. A more complete description including time 
scale fluctuations thus appears necessary. 

That sample to sample fluctuations should play an im- 
portant role both in the statics and dynamics of disor- 
dered glasses is indeed expected from phenomenological 
arguments, e.g. the droplet scenario^i^ which appears 
to describe simpler models such as ^ relatively well, at 
least in low dimensions i2i Let us recall its main conclu- 
sions. In its simplest form, it supposes the existence, at 
each length scale L, of a small number of excitations of 
size ~ L'' above a ground state, drawn from an energy 
distribution of width SE ~ with constant weight near 
SE — 0. While typically the elastic manifold is localized 
near a ground state, disorder averages of static thermal 
fluctuations at a given scale are dominated by rare sam- 
ples/regions with two nearly degenerate minima. For 
example, as a simple but remarkable consequence, the 
(2n)*'^ moment of u fluctuations is expected to behave as 



((«2) _ (y)2)n ^ Cn{T/L')L^^i. 



(7) 



The droplet picture supposes that the long-time equilib- 
rium dynamics is dominated by thermal activation be- 
tween these quasi-degenerate minima controlled by bar- 
riers of typical scale Ul ^ E^ ■ Little is known about 
the distribution of these barriers, but there is some 
evidence22iS2i^ that -0 « 6*. Even a modest distri- 
bution of barriers, however, due to the Arrhenius law 
tl ^ e'^^/-^, yields relaxation time scales with an ex- 
tremely broad distribution as T — + 0. Some probes of 
this broad distribution are the relaxation time moments 
which may be defined in a variety of ways. One begins by 
defining the relaxation time moments in a single sample, 



dt 



(8) 



which, for a particular disorder realization, describe the 
response of the center of mass coordinate to a (spatially) 
uniform force. For n 1 ((t)L) this gives one defini- 
tion of the relaxation time of a single sample. A set 
of disorder-averaged moments may be obtained then by 
directly averaging the above objects, (t")L, giving the 
averaged response of the system: 



dtru. 



q,t 



(9) 



q=l/L 



Alternatively, the distribution (from sample to sample) 
of the unaveraged relaxation time (t) l is described by 
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a second set of moments, In general, one may 

construct many such objects scaling dimensionally as 
t" but with different physical content. Mathematically, 
this is accomplished by averaging arbitrary products of 
the single sample moments, i.e. {V^) l ■ ■ ■ {t'^'^) l, with 
'Yl!j=iPj ~ Ally of these "n-th" moments may be- 
have as ~ ^a(n)UoL /T ^]^-\^ a{n) > n, or grow even 
faster with Inr" ^ , nor is it clear that the differ- 
ent definitions for a given n exhibit the same growth. 
Indeed, we will ultimately find different operators in a 
dynamical field theoretic formulation corresponding to 
each different type of moment, and some indications that 
indeed different growth rates obtain for each of these. 
A theory of these timescales is crucial to understanding 
both equilibrium response and correlations and to near- 
equilibrium phenomena such as creep jiSiii^i^*^ Exten- 
sive calculations are possible in certain zero dimensional 
toy modelsi^ Although some analytical results have been 
obtained for d > 1 within mean-field limits these 
do not include thermal activation over divergent barriers 
(Ul ~ i"). The FRGiSii2*2ai2i on the other hand, ex- 
tended to non-zero temperature j ^^i^^ seems to capture, 
already at the level of the single time scale approxima- 
tion the existence of these growing barriers. However 
until now neither the rare events nor fiuctuating barriers 
have been obtained in this approach. 

In the main part of our study (Sections III and IV) 
we thus investigate how relaxation time distributions ap- 
pear within the FRG. We first show that the equation 
of motion generates under coarse graining a random 
friction term r]{r) (equivalently a random relaxation time 
r ^ 1/77). It is then natural to define, as a toy model, a 
random friction model which, in the absence of pinning 
disorder possesses a manifold of fixed points indexed by 
the full coarse grained probability distribution of the fric- 
tion. A highly non-trivial question is how this distribu- 
tion will flow under RG due to feedback from non-linear 
terms when pinning disorder is reintroduced. We con- 
sider this question at the one loop level, in two stages. 

In Section III we present a highly simplified analysis, 
but with the merit of explicitly exhibiting the broadening 
of the barrier distribution and allowing for some analyt- 
ical expressions. It yields asymptotically a log-normal 
distribution of relaxation times, i.e. a nearly gaussian 
distribution of effective barriers centered around and 
of typical fiuctuations ~ L^^"^. Such a log- normal tail cor- 
responds to the moment exponents a(n) = 2n^ — n. This 
is compared with numerical results22i22i2£ in the case of 
a directed polymer, where the width was fitted to '--^ L^. 
The question of the width is important since a width 
^ L^/^ is not expected to be large enough to modify the 
creep exponent, as can be seen from reexamining the cal- 
culations of Ref. while a ^ widthS£ would pose 
a problem to this order. We derive, within the same 
approximate scheme (Appendix IE|) . closed equations for 
correlations and response functions (the fact that the 
width grows very fast can be exploited in a resumma- 
tion of the fastest growing terms in the dynamical part 



of the Martin Siggia Rose functional). We find that the 
broadening is sufficiently fast to invalidate some of the 
previous analysis, e.g. the existence of a i/r^ scaling 
regime) . 

Section IV contains the full systematic analysis of the 
running dynamical effective action. It is found that re- 
laxation times distributions are determined by a closed 
hierarchy of FRG equations, each level p corresponding to 
an increasing power of frequency to^ can be solved inde- 
pendently of higher ones, the lowest one being the statics 
p = 0. This hierarchy involves functions parameterizing 
the local cross correlations between pinning disorder and 
random relaxation times. The previous approximation 
corresponds to projecting these FRG functions to their 
values at zero, while in fact the full set of non-linear dif- 
ferential equations obeyed by these functions need to be 
solved, a formidably complex task. A thermal boundary 
layer ansatz (TELA) appears as a consistent solution. It 
extends the one discovered in the statics which was shown 
to reproduce droplet theory type behaviour in thermal 
fluctuations. Here it yields a natural growth for mo- 
ments of relaxation times measured by non-trivial expo- 
nents a{n) ^ 2n^ — n determined by eigenvalue problems. 
Although perturbative control remains a challenge, the 
structure of the dynamical TELA which encodes for bar- 
rier distributions opens the way for deeper understanding 
of the field theory approach to glasses. 

The detailed outline of the paper is as follows. In Sec- 
tion |n] we recall the standard results of the FRG for the 
equilibrium dynamics using a single relaxation time ap- 
proach. We then give a qualitative derivation of the two 
scaling regimes for the equilibrium response and correla- 
tion functions. The detailed equations obeyed by these 
functions are derived using a Wilson scheme in Appendix 
IXI and their analytical form is analyzed in the equilib- 
rium regime in Appendix ^ and in the aging regime in 
Appendix In Section UTTl we go beyond the single re- 
laxation time approach. The random friction model is 
introduced in Section flll El We then incorporate pinning 
disorder in an approximate way in Section [HI Gl analyze 
the resulting distribution of relaxation times and show 
that it become broad. The breakdown of the ujTk scaling 
is analyzed in Section IIII Dl In Section IIVI we discuss 
the systematics of the structure of the dynamical effec- 
tive action. It does contain the statics which its recalled, 
together with its thermal boundary layer ansatz solution, 
in its equilibrium dynamics formulation in Section [IVAI 
Then in Section IIV ElIVEl we display the hierarchical 
structure of the FRG equations and how a generalized 
thermal boundary later structure appears as a consistent 
solution determining the growth of the moments of the 
relaxation times through non-trivial eigenvalue problems. 
We conclude in Section^with some general remarks. Fi- 
nally, a set of appendices elaborate on various details and 
calculated related to the main text. 



4 



II. SINGLE TIME-SCALE APPROXIMATION 



correlations are given by 



At conventional pure critical points, scaling emerges 
directly from the existence of a RG fixed point. More- 
over, in an epsilon expansion, the scaling functions are 
obtained to leading order by a simple matching proce- 
dure (REF). The static equilibrium FRG for the ran- 
dom elastic problem is, aside from the complication of a 
functional fixed point, very similar to such an ordinary 
RG calculation. The important distinction is the non- 
analyticity of the T = fixed point function A* (u) which 
at finite temperature results in a narrow boundary layer 



for small 



< 



Tie, whose width continuously decreases 



under the FRG as the running effective temperature T/ 
(see below) flows to zero. The corresponding growth of 
the mean-squared curvature of the effective potential felt 
by the manifold is a hint of unconventional behavior not 
present in ordinary critical theories. 

The influence of this divergence is very dramatic in the 
dynamics. In this section, we attempt to naively extend 
the conventional RG approach to calculating response 
functions to the random manifold problem. This con- 
ventional approach implicitly assumes the existence of a 
single time-scale (at a given wavevector), as we shall soon 
see. This assumption leads to somewhat unsatisfactory 
results for the response function, forcing us to reconsider 
the distribution of time-scales in Section Iflll Although 
we will ultimately conclude that the single time-scale cal- 
culation is fundamentally incorrect, it is useful to review 
the methodology of this approach. 

We begin by reviewing the basics of the FRG. We 
employ the Martin-Siggia-Rose (MSR) formalism,'^'' in 
order to use field-theoretic techniques. The MSR ap- 
proach is based on the generating functional Z[h, h] for 
the disorder-averaged correlation and response functions. 



Z[h, h]= DuD 



(10) 



where the dynamics in ^ is encoded in the action 
S[u, u] = Sq[u, ii] + Sint[u, u], with 

So[u, u] ^ I iuj-t (rjdt - V^) Urt 



-r]T / {iUrt){iUrt) (11) 
Jrt 

Sint[u,u] ^ -\ I {iUrt){iUrt')^{Urt-Urt') (12) 



where h, h are source fields, and we have put a overbar 
on the friction coefficient -q ^rj to indicate the mean, i.e 
that it is for now a constant uniform in space. As justified 
below we have set c — 1. We use the Ito convention to 
regularize equal-time response functions, i.e. Rq{t,t) = 
and Rqit^ , t) = l/rj. The disorder averaged response and 



Rq{t,t') 



S{uqit)) 



Shqit 



(13) 



Cq{t,t') = {Uq{t)u^q[t')) - {Uq[t)u^ q{t' )) S . (14) 

The FRG in its Wilsonian formulation begins by in- 
troducing an ultraviolet (short distance) cut-off A on the 
spatial Fourier wavevectors. In the FRG, this cut-off is 
progressively reduced to A; — Ae~' {Q < I < oo). At 
each stage of the RG, the spatial Fourier components of 
u, u are divided into "slow" and "fast" modes, with mo- 
menta in the range < /c < kie~'^'' and Aje"''' < fc < A;, 
respectively. The fast modes are then integrated out, 
working perturbatively in Sint to one loop order, and I 
is increased by dl. This leads, in the limit of infinitesi- 
mal rescaling dl ^ 0, to a smooth renormalization of the 
effective action for the remaining slow modes, and hence 
to continuous FRG equations for the running disorder 
correlator Ai(u). 

Naive power counting (see e.g. Rcfs.^^'^^) indicates 
that all terms beyond those in 1)1111 - 1)12(1 are irrelevant, 
so we for the moment neglect their generation under the 
FRG. The fiow of the random pinning correlator /^i{u) 
has been derived many times previously .igiigiSSiSS It is 
better expressed in terms of the dimensionless rescaled 
pinning force correlator A;(u) defined such that: 



A^u) = ^e-''e^^'Ki{ue-^') 
Ad 



(15) 



with Ad = 5'rf/(27r)'^ = 1/ {2'^-^TT'^/^V{d/2)), and reads: 

diK{u) = (e - 2C)A(w) -f C"A'(u) + fjA"(u) (16) 
+K" {u){K{{)) ~ K{u)) ~ K' (uf (17) 

Here the fluctuation dissipation theorem ensures that the 
temperature T/ = T is uncorrected but the effective di- 
mensionless temperature T/ = ^dTA''~^e~^' itself flows 
to zero, controlled by the energy fluctuation exponent 
9 = d — 2 + 2(, the temperature being formally irrele- 
vant. Here and in the following we do not not make any 
spatial or temporal rescalings of coordinates or momenta. 

Study of the one loop FRG equation shows that, with 
the proper value for the roughness exponent C ~ 0(e) de- 
pending on the universality class (RB,RF or RP), the di- 
mensionless disorder correlator converges non-uniformly 
to a nonanalytic "fixed-point" function A*(u) formally 
of order ~ 0{e) as / oo, whose functional form is 
not important for this discussion (see however Section IV 
for much more details). The non- uniformity of this con- 
vergence is due to a boundary-layer centered on m = 0, 
whose width decreases continuously with scalem^^ In 
particular one can show that, to this orderiiSiS^: 



hm r,A;'(o)^-x' 



(18) 



Thus asymptotically the curvature of the correlator di- 
verges with the scale (here x = |A'*(0^)l is a constant 
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depending of the universality class, e.g. for periodic sys- 
tems X = e = 4 — d. 

To one loop, all that remains is the renormalization of 
the mean friction coefficient, since the elastic modulus 
is fixed by Galilean invarianceiiiSiS^ and the temper- 
ature by the Fluctuation-Dissipation-Theorem (FDT). 
This was determined in Refs. 



(19) 



where T/ 
scale as: 



-Ar(o) 



■^i-^+oo X^/Ti thus grows with the 



(3 = T*/T 



(20) 
(21) 



where (3 is the reduced bare inverse temperature and 
T* = x^A?~'^ /Ad a non-universal temperature scale. fH^ 
implies activated scalin^^t^ since the friction coefficient, 
which plays the role of a time scale r = rj/ h?, grows ex- 
ponentially with the length e': 



m = exp 



(3 



(22) 



(in d = 4 one has Ti = pe'^^/P and /3 = x^A-V(TA4)). 

Activated dynamics leads to ambiguities in a single 
time-scale approach, as can be seen from a simple match- 
ing argument. We consider for simplicity the equilib- 
rium dynamics, in which the response and correlation 
functions are time-translationally invariant (TTI). It is 
then convenient to work in terms of both frequency lj 
and wavevector. The usual RG considerations lead one 
naively to expect that the response function obeys a re- 
lation, 



i?fc(w) = e'^'RkeiiujTi), 



(23) 



where ti = e'^^rji/rjQ. We now obtain two inequivalent 
scaling forms by matching. In particular, if we choose 
fee' = 1 (we set A = 1 for now), we find 



1 



(24) 



X and Y scaling limits. The first is formally defined by 
defining the scaling variable Y = ujTk . With Y fixed and 
cj,/? ^ 0, one obtains the scaling regime in H24|l . The 
second scaling regime obtains with X — fc^|lna;| fixed 
and cij, fc — > 0. To reconcile the two regimes, note that 
\nY - {x/9 + X)/k\ so that for fixed X, as cj,fc 
(the X scaling limit) Y ^ oo\ The second scaling regime 
((|2SJ|) thus appears to occur at the "boundary" {Y = oo) 
of the first. 

We have indeed verified this behaviour, and obtained 
the analytical scaling forms similar to those in (|24I2()|) 
using FRG techniques under the assumption of a single 
characteristic time scale parameterized by rj. These cal- 
culations are performed in Appendix A for equilibrium, 
and in Appendix B for the more general noncquilibrium 
situation. Interestingly, the general equations for the re- 
sponse function in this approach share some similarities 
to those arising in the (infinitely connected and/or large 
N) mean field limit of a number of model glasses. As 
discussed in Appendices A and B their solutions exhibit 
several time regimes with various aging scaling forms and 
also show differences compared to mean field. 

Many features of these results, however, point to prob- 
lems with the single time scale assumption, as also dis- 
cussed in Appendix A. The real-time response function 
is found to be an increasing function of wavevector at 
fixed time in the logarithmic {X) scaling regime. This 
somewhat unexpected (and possibly unphysical) behav- 
ior is apparently a very general consequence of the mere 
existence of two distinct scaling limits, and hence is in- 
evitable given the single time scale approach. More sig- 
nificantly, the appearance of a sharply-defined in the 
mean response function (in the Y regime) is difficult to 
understand on physical grounds. Even in (random) mod- 
els involving only a small number of degrees of freedom, 
while a given sample may be characterized by a longest 
relaxation time, the sample to sample variations of this 
would generally lead, as espoused in the introduction, to 
the disappearance of such a time in the mean response. In 
the collective elastic model considered here, interactions 
between the enormous number of modes with differing 
wavevector (and hence differing relaxation rates) would 
only worsen the situation. 



with 



1 ^(k~ 



(25) 



If, however, we choose wt; = 1, we find asymptotically 
i?f = (iln(i))2/«7^(2)(ifc^|ln^|). (26) 

Note that these two forms cannot be related by redefining 
the two scaling functions VS^^'^^ , as can be seen from the 
fact that In(wTfe) ~ \nw+[x/0)k-'^ = {x/ 6+k'^\\nuj\) /k\ 
which is not a function of Inwj only. 

The inequivalence of H24I26|I appears to point to the 
existence of two scaling regimes, which we will call the 



III. BROAD DISTRIBUTIONS OF TIME 
SCALES: SIMPLIFIED APPROACH 

A. distribution of relaxation times and the f-term 

Up to this point, we have assumed that the dynamics 
at each scale can be described by a single friction coef- 
ficient Ti = rji, which corresponds to a sharply defined 
time scale for relaxation. On general grounds, however, 
we should expect extremely broad distributions of relax- 
ation times. This follows simply from the Arrhenius law. 



Tk = T;=in(fc/A) roexp(C/i;/r), 



(27) 
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which estimates the time required to overcome an energy 
barrier of height C/fe at scale k. At low temperature, even 
a modestly wide distribution of C/fe gives rise to extremely 
broadly distributed r^. If this distribution is sufRciently 
broad, it is no longer adequately characterized by its av- 
erage, and indeed many physical quantities may depend 
upon the precise form of the distribution. 

We now investigate how this distribution can be in- 
corporated into the FRG treatment within the MSR for- 
malism. We will consider a spatially varying friction co- 
efficient 77(r), with the equation of motion modified 
to 



r]{r)dtUrt = cVlurt + f{urt,r) + C{r,t), 



(28) 



where, in order to maintain the stationary equilibrium 
Boltzmann probability distribution function, the noise 
correlations are modified to 

(C(r, t)C(r', t')) = 2ij{r)TS''{r - r')6{t - t'). (29) 

For simplicity, we will initially take r]{r) to be identi- 
cally and independently distributed at each r, according 
to the distribution P{r]). The distribution naturally en- 
ters the MSR theory via its characteristic function, which 
we parameterize by F{z): 



+ 00 



(30) 



The Taylor series expansion of F{z) thereby gives the 
connected cumulants of rj, 



oc 

E 

m— 1 



/ 1 \m,+ l 



F\z\ ^ > rj'-'"'- z" 



(31) 



where ry*^™^ = [?7™]c is the m"^ cumulant (connected) 
moment of rj. For the continuum field rj{r), the analogous 
expression is 



exp{— I 7j{r)z{r)) = exp{— / F[z{r)]) 



(32) 



We initially assume no cross-correlations between 77 (r) 
and f{u,r), though these can to some extent be gener- 
ated in perturbation theory. The single time scale model 
studied in the previous Section (and Appendices A and 
B) with ri{r) = rj corresponds to F(z) — rjz. Although 
this is not essential, we shall assume here an initially 
narrow (but not (5- function) distribution P{ri). As shown 
in Appendix C, even if initially F = rjz, a higher-order 
analysis shows that a non-trivial distribution is generated 
under coarse graining. 

In the MSR formalism, the modified equation of mo- 
tion, (|28I29I) . is described by the action: 

So[u,u] — / rj iUrtdtUrt — iUrtC^'lurt — rjTiUr-tiUj.t, 
Jrt 

(33) 



Sint[u,u] = J F[J {iUrtdtUrt - TiilrtiUrt)] 

{iUrt){iUrt')^{Urt - Urt')- (34) 



We have defined F[z] — rjz + F[z] where F[z] starts 
with higher powers of z, to do perturbation theory using 
the average friction coefficient. One immediate remark 
is that the statistical tilt/translational symmetry (STS) 
holds, thus c will not be corrected and so we set it to 
c = 1. Note that the rj kinetic term could equally well 
be considered as an interaction term, in the spirit of a 
"perturbation theory in iw" with bare propagator sim- 
ply = (in real time Rq{t,t') = l/q'^S{t- - t')). 
Note also that for each realization, the instantaneous re- 
sponse function satisfies 



R{v,v',t-t' = 0- 



77(r) 



5{y - r'). 



(35) 



Averaging over disorder gives Rk{t — t' = 0+) = l/rj. 

Although it may seem obvious, it is important to stress 
at this stage that the renormalized relaxation time mo- 
ments are measurable quantities. One can define the 
renormalized moments in the usual way from the effective 
action F. Taking the same form as 134(1 . one has 



(2) 



UrtUrtUrt'Urt 



(36) 



with of course many higher order terms describing the 
higher friction coefficient moments, momentum depen- 
dence of vertices, etc. As usual in field theory, correla- 
tion functions are exactly evaluated at tree level using 
this effective action - thus the 77^^^ term here has the 
meaning of a fully-renormalized second moment on the 
scale of the system size (or infrared momentum cutoff). 
One may then consider the physically defined relaxation 
time from (jSl and construct its second moment: 



(*)i - L2{d+2} 



dtdt' tt' I {UritUr'^o){Ur2t'Ur'^o)- 

(37) 



On physical grounds, in equilibrium, we expect that the 
latter product of response functions in two "replicas" is 
the same as considering the product of two subsequent 
responses in a single replica, provided the two response 
measurements are taken far apart. 



= lim (Urit+TUr'^TUr2t'Ur'^o)- 



(38) 



This latter four-point function can be calculated using 
the effective action above. One finds 



lim 

r — ^00 



{Ur^t+rUr[TUr2t'Ur'^a) (39) 

rir[r'^r2 
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FIG. 1: Top: graphical representation of the disorder vertex 
(double lines) and of the F term vertex (the dots represent 
the time derivatives) where the arrows represent a u response 
field and solid line a u field (arrows are along inreasing time) . 
Bottom: (a) corrections of the F term to itself at T = which 
vanish (b) corrections to order F^ which also vanish (as all 
orders do - see text). 

where go ~ 1/^ is the infrared momentum cutoff, and * 
denotes the convolution in the time domain. Integrating 
over the time coordinates and using the result of statis- 
tical translational symmetry i?(qo,0) = 1/Qoi one then 
obtains 

WL-ivL^f+V^'^L^-"- (40) 

For a not too broad distribution of friction coefficients 
with scale-independent rj^'^\ the correction due to 77^^^ 
is vanishing for c? > 4, and small compared to the first 

"disconnected" term (scaling as {t)L ) for any d. How- 
ever, for the glassy dynamics studied here, we will find 
77^^) is exponentially larger than if' as a function of L, so 
that in fact the second term is dominant. T hus t he sec- 
ond moment of the physical relaxation time {t)\ indeed 
measures the coupling constant rj^'^^ as promised. 

B. no pinning disorder: the random friction model 

We now turn to the FRG analysis of the modified ac- 
tion in H33I34|I . We first consider only the effects of ran- 
domness in rj, neglecting the pinning disorder A. This 
defines a random friction model described by the MSR 
action with A = 0. Remarkably, the random friction 
model represents an infinite manifold of fixed points pa- 
rameterized by F[z\. Indeed, a diagrammatic treatment 
explicitly shows the absence of renormalization of F[z\ 
order by order in F . Despite this absence of renormaliza- 
tion, the random friction model represents a non-trivial 
interacting field theory. 

For simplicity, we sketch this here for T = 0. In this 
limit, the vertex is 

J^F[JjUrtdtUrt]. (41) 

Diagrams occuring in the expansion of F are indicated 
in Fig. n (a) and (b). The fields connected by dotted 



lines occur at the same spatial point, solid lines with 
and without arrows indicate u and u fields, respectively. 
Considering a product of the form 

/ (/ «Writil9tiUriti)"' / (/ mr2t2^t2Wr2t2)"^ (42) 

the only possible contractions contain products of the 
type < dtjUrjtjiUrktk > with no time loop allowed. Thus 
all relative time integrals factor and one is left with prod- 
ucts of integrals of the type dtR{rj — rk,t) which van- 
ish since the response function vanishes for i < and 
t — > -l-oo Thus F does not correct F. However F 
itself produces new terms such as: 

/ iUrtdtUrt (43) 
Jr,t 

I iUrtidf^UrtiiUrt2dt2Urt2, (44) 

r,ti,t2 

obtained by time gradient expansions, with nonvanishing 
coefficients (of the form ^ tdtR{rj — rk,t)), as well as 
similar terms with higher order time derivatives (note 
that similar terms containing also higher order spatial 
gradients are also generated, but we will not consider 
them as important here^). One can embed these new 
terms into a new function: 

J^F2[JjUrtd^Urt] (45) 

and so on - the full systematics of these new terms will 
be examined later in Appendix |E| and Section (jIVp . It 
is important to note for consistency that there is also no 
feedback from higher-derivative terms such as F2 back 
into F. Graphically, as in Fig. (a) and (b), one can 
perform time integration by parts along each line joining 
several vertices which leads to terms dfutdfUt with p+ 
q > 1. Indeed a very useful rule is represented graphically 
on Fig 121 one can simply shift the time derivative along 
any internal response line to the external one (at T = 
any diagram is a tree of such lines) since, schematically: 

Ut{iutdtiUtj^)iuti ^ utdt^Rt^^tiuti (46) 
= -utdtRt^tiuti dtUtRtutiuti (47) 

after integration by parts. In the Fourier domain, this 
rule is just conservation of frequency along all solid lines, 
since the interactions are all fully non-local in time and 
therefore do not carry frequency. 

To conclude, the apparent non renormalization of F[z\ 
makes it tempting to define a manifold of fixed point 
theories indexed by F[z\. These fixed points are quite 
interesting and non-trivial. For instance, the computa- 
tion of the averaged response function at T = can be 
mapped exactly onto the problem of calculating the par- 
tition function for a self-avoiding walk. This is developed 
further in Appendix D. 




FIG. 2: Shifting of a time derivative from an internal line (in 
the middle) to an external one (right) along a line of response 
functions at T = (works for disorder as well as F vertices) 

C. pinning disorder: distribution of barriers 

We now consider the combined effects of the pin- 
ning disorder and distribution of time scales. Because 
the pinning disorder can be defined in a purely static 
theory (using the equilibrium Bolzmann partition func- 
tion), its renormalization is unaffected by the F term. 
However, the converse is not correct. Due to the non- 
renormalization of F in the random friction model, we 
must consider only terms of order F^A'^, with q > I. 
The leading non-vanishing terms correcting F at 0(A) 
are linear in F, and are indicated diagrammatically in 
Figs. I3I4I They are computed in detail in Appendix El 
but one easily sees the structure of the result, thanks to 
the property of shifting internal time derivatives (dots 
in the figures) to the external ones, e.g. that the three 
graphs in Fig. |3| have identical values. 

There is a subtle distinction between the contributions 
in Fig. O and those in Fig. 0] In particular, in the dia- 
gram of Fig. 01 the pinning vertex suffers contractions be- 
tween both of its independent time variables (i.e. graph- 
ically both ends of the double-line) and the same time 
variable (dotted leg) of the f-term vertex. The locality of 
the response function therefore implies that the two inter- 
nal times of the pinning force correlator are constrained 
to be nearby, justifying a temporal gradient expansion 
and hence giving a leading contribution proportional to 
A"(0). In the diagrams of Fig. O by contrast, the two 
times of the pinning vertex are contracted with different 
legs of the F-term. These diagrams therefore generate in 
fact more general terms involving A" (ut — Uf ) with free 
integration over t and f. If \ut — Ut'\ is not extremely 
small (within the boundary layer), this is a small correc- 
tion lacking the singular temperature dependence. It is 
thus not clear at this stage whether or not the graphs in 
Fig. should in fact be interpreted as renormalizations 
of the F-term. 

In fact, the true situation is more delicate, and will be 
returned to in Sec. (|IV|l . For the moment, however, we 
will shut our eyes to this complication, and gain some 
physical insight by taking into account both sets of di- 
agrams as renormalizations of the f-term. Their sum, 
integrated in the momentum shell, gives the following 
correction to F: 

SF[z] = ~A"{0)SdAf-UlizF'[z] + 2z^F"[z]). (48) 

Here the F' and F" terms comes from the diagrams in 
Fig. I4I3I respectively. In agreement with H19I) it can be 
rewritten as 

diFi[z] = (9, lnr],)izF;[z] + 2z'F!'[z]). (49) 
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FIG. 3: Graphs involving pinning disorder which correct the 
F term proportionally to F" 




FIG. 4: Graphs involving pinning disorder which correct the 
F term proportionally to F' 

Note that the disputed diagrams in Fig. O do not con- 
tribute to the mean relaxation time due to the second 
derivative of z, so that rj is unambiguous. The RG equa- 
tions for the connected moments, ri^'^\ of 77 are thereby 
obtained using as 

vl"^ -vi"\^ f (50) 
Vo 

It is more convenient and physical to introduce the 
random barrier U — Ini], and the barrier corresponding 
to the average relaxation time Ui = Inr^j. Changing to 
the energy variable v = Inz, and letting Gi{v) = Fi[e") 
gives 

dufii[v]^2G'l[v]~G[[vl (51) 

i.e a diffusion with drift equation. Some physical under- 
standing of G{v) can be obtained from the two extreme 
limits, 

Gw-{r-inp(o) ' (52) 

as can easily be found from assuming a constant 

probability density for small barriers, < P(0) < 00. 
More generally, using (|30|l the diffusing and drifting "den- 
sity" Gi is related to the barrier probability distribution 
via 

Gi[v] = - In / dUPi{U)e-''"^\ (53) 



9 



Formally, the solution of H51|l is given by 



Gi{v) = dw 



exp 



jv-Ui- w)^ 
8Ui 



(54) 

and inverting the results via H53|) to obtain Pi{U). A 
simple approximation may be applied in the regime of 
large U Ui and [/; > 1, in which Gi{v) < 1. In 
this case, it is valid (and justified a posteriori since the 
distribution of barriers become broad) to replace e"*^ 
by e{U < -v) and thus one gets that Pi{U) « GJ(-[/). 
This yields (via H54|l or directly differentiating (ISlfl l 



PiiU) 



1 



exp(- 



8Ui 



), U>Ui. (55) 



Note that this asymptotic form reproduces all cumulants, 
rji"^ ~ (exp(nC/))c ~ exp((2n^ — n)Ui) as expected. 

H55|) is clearly not exact. Indeed, a breakdown of H55|) 
is inevitable on physical grounds, since the mean/typical 
barrier cannot be negative! It suggests a distribution of 
barriers with a width proportional to \/Ui, and hence 
a peaked distribution (since the mean barrier oc Ui). 
Nevertheless, it does represent a very broad (in fact log- 
normal) renormalized distribution of characteristic times 
ij. While there is no reason to believe that such a log- 
normal tail is exact, the true distribution of relaxation 
times will certainly be very broad, with significant con- 
sequences for the average response functions. 



D. Breakdown of uTk scaling 

The first consequences of this broad distribution occur 
in the variance 77^^^ of the relaxation time, and hence at 
0{uj'^) in the response function. We therefore examine 
more carefully the O(w^) terms in the dynamical action, 
but for the moment still neglecting the full functional 
dependence of these terms (i.e. on ut — ut'). In the 
kinetic part of the action (representing relaxation times 
and their fluctuations), we include the following terms: 

^kin ^ [V iUrtdtUrt + D iilrtdtUrt] 

Jrt 

-trt2- (56) 



77(2) f 



One has in general, defining (i")^ — t"-Rk{t)/ Rkit)'- 



2D 



(57) 
(58) 



There is no generic constraint on the sign of D. If the 
inverse response function contained only the two above 
terms {rj and D), then causality requires D to be positive 
( similar to an inertial term) .^"^ Since in general these are 
only trucation of an infinite series of terms in power of 
iw, the only constraint is causality, i.e. that all poles 



FIG. 5: Correction to the (icj)^ term in the response function 
coming from the second moment r/'^-* of the relaxation time 
distribution 



in uj lie on the same side of the real axis. These three 
couplings satisfy the following closed RG flow equations 
to first order in A: 



diDi 



(2) 



(2) 



6F/?7, 



(2) 



(59) 
(60) 

(61) 



where F; = — A"(0) ~ /3e^' and the correction to D from 



(2) 



is the graph represented in Fig. (0). 



77 

For 77*^2) = the equations for Di and 77; , which can be 
obtained e.g. from (|A2IIA6(I by expansion to second order 
in iuj, are consistently solved with Di ^ A^^rjf, in the 
limit of large I, where rji = rj^ exp[/3(e^' — 1) /6], consistent 
with the Taylor expansion of the putative scaling function 
g{y = iujTk) given in (jA9|) based on the single time scale 
analysis (|A9I) . 

For 77(2) > 0, however, the broad distribution of re- 
laxation times completely alters the situation. From 



the above equations, 

J2) ^ -2 



(2) 



^70 



and hence 



rjl ' ^ 77f (the mere exponential prefactors are negligi- 



ble) at large I. Thus the feedback of 77^2) d dominates 
the renormalization of Z?, and at large / one finds: 



d-1 



77o 



(62) 



Thus, allowing for fluctuations in relaxation times in- 
validates the WTfe scaling form already at order ijp'X 

It is still possible, within the approximation scheme of 
the present Section, to obtain an equation for the disorder 
averaged response function. This is explored further in 
Appendix |E| 



IV. DISTRIBUTION OF TIME SCALES: FULL 
STRUCTURE OF THE DYNAMICAL FIELD 
THEORY 

We have established the mechanism for breakdown of 
the unphysical lot-^ scaling regime and described the in- 
dications of a broad distribution of timescales within the 
FRG. However, to properly determine this distribution 
and its consequences, e.g. on the mean response func- 
tion, requires a much more complete analysis. While we 
have unfortunately so far been unable to carry this pro- 
gram to completion, in this section we will detail the for- 
mal structure within which this analysis must be carried 
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out. In particular, we shall see that the distribution of 
relaxation times and its consequences is encoded within 
the boundary layer (BL) regime already present in the 
statics. An understanding of the equilibrium dynamics 
is therefore contingent first upon an understanding of the 
static BL, and we first describe the rather complex struc- 
ture therein. Following this discussion, we show how the 
BL regime recurs in dynamical theory, and show how 
it can be formulated to describe broad distributions and 
non-trivial scaling of the moments of the relaxation time. 



tity. We can therefore benefit from the knowledge of the 
thermal boundary layer in the statics. To do so let us 
review how the static disorder correlations are encoded 
in the dynamical formalism. 



A. statics thermal boundary layer 

In the appropriate limit the dynamical theory should 
reproduce the results for the corresponding statics quan- 



At equilibrium, the part of the dynamical action con- 
taining the static disorder correlations comprises those 
terms with no explicit time derivatives, and reads: 



Smt=-T: f iUrtiiUrt2^{Urti - Urt2) - 7: [ iUrtiiUrt2'i'U'rt3Sli^\urti , Urt2 , Urt^) - ■ ■ (63) 

This form is easily understood as arising from the cumulants of the pinning force. The relation was given for the 
second cumulant in (j^J and for higher ones it reads: 

/K,ri)../K,r,f = {-fS^^\u^, ■ ■ ■.uk)5\r^, • • • , r^) (64) 

with sf^ {u, u') = A{u—u') as Q. Dues to statistical translational invariance S'^'^^ (wi, • • • , Uk) — S^'^'' (ui+v, • • • , Uk+v) 
and satisfy reflection symmetry S'--^\—ui, • • • , —Uk) = (— )'^S'*-'^-*(wi, • • • ,Uk)- The cumulants higher than second are 
generated by coarse graining, and are thus included here from the start. 

The static problem being defined from the equilibrium Boltzmann measure (cflSJ, deals not with the distribution 
of the random force but with that of the random potential 

V{ui,r,)..V{uk,rkf = S^'^^u^, • • • , Uk)S''{n, rk). (65) 
Since f{u,r) = —duV{u,r) one has: 

A(u) = ~R"{u) (66) 

S^f\ui,U2,U3) ^ did2d3S^^\ui,U2,U3) (67) 

and so on. 

It is straighforward to derive the one loop FRG equation in the Wilson scheme for these cumulants using the 
dynamical formulation. They are conveniently expressed using rescaled cumulants: 

and read, for the second and third cumulant: 

diMu) = (e - 2C + Cudu)A{u) + fiK"{u) + 2^ioo(0, u, 0) - K\uf - K"{u){K{u) - A(0)) (68) 
diS{ui,U2,U3) = (-2 + 2e - 3C + Cui9„JS'(ui, U2, U3) + -T;(S'2oo(wi, "2, "3) + 'So2o(wi, -"2, "s) 

+ <S002(W1, U2, U3)) - ^A(0)(S'20o("l, U2, U3) + ^02o(mi, M2, U3) + 5o02 (^1 , U2 , M3) ) 
-^A{ui - U2)Siio{ui,U2,U3) - ^A"(mi - U2){S{ui,Ui,U3) - S{ui,U2,U3)) 

-^A'{ui - U2){Soio{ui,U2,U3) - 5ioo(wi,M2,W3) + Sqio{ui,Ui,U3) + Sioo{ui,Ui,U3)) (69) 



I 

where we have denoted S^^P — S and we have suppressed explicitly the feedback of the fourth cumulant S*^**' into 
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the third one. One can check that this gives exactly the 
derivatives IjSTI ) of the one loop FRG equations for the 
static correlators R and S'^'^^ displayed in (6,7) in Reff^l. 
These relations (|67l ) should indeed be preserved by RG 
at equilibrium. 

As discussed in Rcf. 23 when all arguments of these 
functions are distinct and order one conventional scaling 
holds. That is, at large scales for which T; ^ 0, the func- 
tions 5^'^-' approach well defined nonanalytic fixed point 
forms 5''^^*. Moreover, naively these can be organized 
in an e = 4 — d expansion in which 5'('=)* ~ e'^, fc > 3. 
Naively this would allow the truncation of the hierarchy 
of FRG equations for the S''^\ neglecting feedback of the 
k > p cumulants with an accuracy of O(e^). However, 
the convergence to these values is highly non-uniform as 
mentionned in Section II since at non-zero temperature 
these functions remain analytic at w = 0. A detailed 
analysis of the static hierarchy of FRG equations relat- 



where Sd are well defined functions of order one in the 
TBL u ^ 1. The set of (^-dependent) constants, 

/2, = ^^'^^(0,---,0)/(xe)^P (73) 

with /2 = A(0)/(xe)'^, have the meaning of the lin- 
earized random force cumulants within the zero temper- 
ature Larkin description. As discussed in Ref. 23 the 
crucial difference with the naive dimensional reduction 
result, where the /2p are unrenormalized, is that they get 
feedback from the TBL functions and acquire non-trivial 
asymptotic values. 

The TBLA encodes a huge amount of physics - in par- 
ticular, all the distributions of minima degeneracy re- 
sponsible for large thermal fluctuations in the droplet 
picture, as detailed in Ref. I^^- For instance, averages 
such as (UI) can be estimated using the TBLA, the coeffi- 
cients c„ being in principle determined by the functions 

This already non-trivial structure must now be gener- 
alized to intrinsically dynamical quantities. 

B. dynamical hierarchy of kinetic coefficients 

In a conventional dynamical renormalization group in 
the MSR formalism a succession of individual terms are 
added to the action corresponding to increasingly high 
frequency kinetic coefficients, e.g. for a particle the 
Stokes drag, inertial mass, ... For the disordered elastic 
manifold however we recognize that these kinetic coeffi- 



ing these cumulants revealed the existence of a thermal 
boundary layer (TBL) of the form: 

Aiu)^A*iO)-fifiu) (70) 
u = exu/fi (71) 

for u — 0(1), Ti <^ and / an analytic function with 
f{x) ^ \x\ at large x to match the cusp of the zero tem- 
perature solution. For higher cumulants the very uncon- 
ventional TBL scaling implies that it is no longer legit- 
imate to neglect the feedback of higher cumulants (the 
n-th cumulants gets a feedback from the n and n + 1 
ones). Therefore, we are unable to truncate and solve 
the hierarchy of FRG equations. Instead, in Ref. we 
argued for the consistency of a thermal boundary layer 
ansatz (TBLA), which for the force cumulants reads: 



(72) 

I 

cients have a broad distribution characterized by an infi- 
nite set of cumulants and cross-correlations, which more- 
over can be non-trivial functions of displacement field 
differences. The latter dependence was neglected in the 
approximate treatment of Section III. The need for treat- 
ing it was already indicated in the ambiguities in the 
diagrams of Fig. 3. Each of these cumulants and cross- 
correlations appears as a distinct interaction function in 
the MSR action. 

By symmetry (time translation and STS, statistical 
reflection, causality) alone, the set of all such interactions 
contributing to the effective action at zero temperature 
can be written as 

oo „ 

4-cxD n 

S^^\urt„..UrtJ n Uidl^ru f (74) 

fe=l i=l 

where > and from STS symmetry, the S{ui + 
u, ..Un + u) = S{ui, ..Un) are translationnally invariant, 
and statistical reflection implies the full action is also 
invariant under {u,u) — > (— u, — u). The random force 
correlators correspond to 

'5'p"io('"i' ■ • • : = '5'rf"^("i, • ■ ■ , Un), (75) 

where P ~ above indicates the function with p'^ — 
for all i, k. Other terms correspond to intrinsically 
dynamical cumulants. 



5(fe)/„ „ ^ j/fe + (xe) 7^(4 ("!'•• -^fc) ^ even 
I (xe)" ^Tis\'{ui,---Uk) fc odd 
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It is instructive to begin the characterization of such 
terms at T = by Usting all possible forms in order 
of increasing number m of time derivatives and order of 
cumulant (i.e. the number of independent times which 
equals the number of m fields at T = 0). Each term in ff^ 



can be assigned m = X!"=i ^'k=i ^Pi- ^'^^ organisational 
purposes it is convenient to rewrite the action of H74I) in 
a schematic (but transparent) notation, first expanding 
in number of cumulants: 



5 = zui(fc2 + S(ai))ui + 5„t (76) 
J:{s) ^T]s + Ds'^ + .. (77) 

Sint = -^WlW25i2 - ^iUliU2iU3Si23 - ■■ (78) 

2 b 

with s — iu!. Here the subscripts 1, 2, ... refer to different times being independently integrated over in the action at 
same space point r (further integrated on). The S12... are then functions of the ui, U2, ■■■ and their time derivatives. 
We then expand each of these in increasing number m of time derivatives: 

S12 - A(ui2) + {ill - ^2)0(^12) + UiU2A{ui2) + {ill + ul)B{ui2) + (Ul - U2)C(ui2) + . . . (79) 
^123 = S{ui,U2,U3) + -{uiH{ui;U2,U3) +U2H{u2\U3,Ui) +U3H{U3;UI,U2)) +UlU2W{ui,U2;U3) + .. (80) 

(81) 

As discussed above, each new term in (|79I81|) corresponds to statistical properties of the random kinetic coefficients 
and forces in a renormalized equation of motion, in particular, 

h D{u, r)u + r/(M, r)u = V^u + f{u, r) + g{u, r)u'^ H h Q{r, t), (82) 

with 



D{u,r)=D, ij{u,r)f{u',r') = -G{u - u')S{r - r'), r]{u,r)r]{u' ,r') ^ A{u - u')5{r - r'), (83) 

g{u, r)f{u', r'f = B{u - u)5{r - /), /(w, r)D{u\ r'f = C{u - u)S{r - r'), (84) 
c 1 

V{ui,ri)f{u2,r2)fiu3,r3) = -i?(wi; -U2, U3)(5(ri - r2)(5(r2 - rg) (85) 

c 1 

ri{ui,ri)ri{u2,r2)f{u3,r3) = -M^(mi, M2; ""a) (86) 



In the approximate treatment of Sec. lIIII Ty*^^-* hence corre- 
sponds to A{u) approximated as ^(0). Note that it is the 
small argument behavior of A{u) (and its higher cumu- 
lant analogues) that is related to the physically interest- 
ing second (higher) relaxation time moment 77*^^^ ~ (t)'^ 
(7y(") ~ (tj'^ ). Hence these relaxation times are encoded 
within the BL regime of these functions. This was also 
apparent from the naive renormalization of rj by A"(0), 
also a BL quantity. We will return to the problem of the 
dynamic BLs in G{u), A{u), . . . momentarily. 

Although it is convenient as above for the purpose of 
enumerating terms in the dynamical action to first sep- 
arate by cumulant index n and then by number of time 
derivatives m, conceptually we analyze them in the op- 
posite scheme, i.e. collecting all terms of a given m, and 
organising these afterward in order of n. This scheme is 
clearly convenient insofar as the first term (m = 0) of 
each of the Si2---n corresponds to the n-th term of the 



static cumulant hierarchy, so that the set of terms with 
m = satisfies a closed hierarchy of FRG equations in- 
dependent of those with m > 0. We now demonstrate 
diagrammatically that, at zero temperature, a similar 
property holds for m > 0. In particular, all terms of 
any given m will satisfy a closed hierarchy of FRG equa- 
tions containing only terms with m' < m. Thus, one 
may imagine (dream of?) solving the FRG equations up 
to level m, then using this solution to complete a closed 
set of FRG equations for level m+1, and iteratively solv- 
ing for higher and higher m. 

This closure relies on the rule of conservation of pow- 
ers of frequency, established in Sec. IIIII Recall that this 
occurs because at T = 0, the correlation function van- 
ishes, and all contractions take the form of causal re- 
sponse functions. Thus no closed time loops can ap- 
pear in any diagram. This implies that internal time 
derivatives which appear in any diagram appear as fac- 
tors of frequency of some external leg to which they are 
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FIG. 6: Compact notation for a generic vertex at T = 0. 




\ □ 



FIG. 7: One loop diagrams which correct the effective action 
at T = 0: the internal lines contain the full response function 
and the graphs are IP irreducible. Graph (a) is a "tadpole". 
Graphs (b) and (c) (and higher orders) correct terms with 
n > 1, 2, 3 respectively. 



connected. In any case, this rule implies that, since all 
terms in the action have m > 0, terms with m' > m can 
never reduce their number of time derivatives by con- 
traction with another vertex at T = 0, and hence cannot 
renormalizc m- vertices. This is true to all orders for dia- 
grams with any number of loops. The frequency conser- 
vation rule implies in fact a more detailed result. If the 
quadratic terms in S(s) (77, D, . . .) are regarded them- 
selves as coupling constanta^^ (with m = 1,2,...), then 
each term in the FRG equation for any quantity at level 
m is a product of factors for which the total frequency 
level (i.e. ^ rrii for all terms i in the product) is exactly 
m. Thus a static quantity (e.g. A) can renormalize a dy- 
namic one (e.g. G,r]) only in combination (i.e. multiplied 
by) another dynamic quantity, and so on. 

One can also establish a set of rules to understand 
how cumulants with different n are connected in the 
FRG equations. At T = 0, this process is highly con- 
strained, since each contraction involves one response 
function, which removes one u, it is straightforward to 
count the possible connections. We will restrict our at- 
tention to one loop diagrams, anticipating future nonper- 
turbative exploration using the exact RG^SSi^ in which 
only these appear (and in any case only these are con- 
sistently treated in the Wilsonian scheme of this paper). 
The counting is illustrated for such one loop diagrams 
in Figs. I6I7I One readily sees that when N vertices are 
combined in this manner, the resulting vertex which is 
renormalized in the effective action contains a total num- 
ber of independent times (or u factors) n = X]i=i ~ 
due to the N response functions appearing in the loop. 

With these rules in mind, we can describe the structure 
of the FRG hierarchy as far as the feeding of terms of a 
given m,n into other m',n'. We note symbolically by 
S"^ the terms with n response fields and m time deriva- 
tives. The term ^ is Si, the response function is the 
quadratic part of (we note = quad(S'^)) and 
the cumulants 77'^"-' are included in 5". From the above 
discussion, neglecting rescaling terms, the structure of 
the FRG equations reads 



m n+1 

XQ{n) _ Q(n+1) , n{n') Q(n+2-n') q("+3-"'-"") , /'07^ 

m' —0 n' —1 m' +m" <ni n'+'n/'<n+2 



It is straightforward to see that this series contains a 
finite number of terms for any given m, n. Let us suppose 
that an A''- loop term exists, such that each vertex making 
it up has Hi time integrations. Suppose of these N ver- 
tices. Hi > 1 for N' of them and rii = \ for the remaining 

N-N'. Then n = J2f^i n^-N = N-N'+Y.^^i ^t-N = 

X]i=i('^i — 1). Hence at most N' < n. Now the remain- 



ing N — N' > N — n vertices have only one time integra- 
tion. Since there are no allowed local terms without time 
derivatives, these must each have rrii > I, i.e. m > N — n. 
Turning this around, iV < m -|- n, so that the series of 
one loop diagrams terminates at at most {m + nY^ order. 
Clearly from (|87|l . each order contains a finite number of 
terms, so that the one loop FRG equations are finite. 
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1. terms proportional to frequency m = 1 equation for terms with n < 2 for which we need terms 

up to n = 3: 

We wiU now examine level m = 1 and m = 2 of the 
hierarchy. For m = 1 we will restrict to study the FRG 



SI 



(1-3) 



_ 1 /■ 1 

■qiUrtilrt-T: I iUrtiiUrt^iUrti - Urt2)G{Urti - Urt^) - 
rt ^ Jrtxti O Jrtit2t-3 



(88) 



where G{—u) = —G{u). 

The renormalisation of rj and G is determined by a standard if cumbersome one loop calculation performed in the 
Appendix^ The corresponding graphs are represented in Fig. |Sland|^ From dimensional analysis and the structural 
form of the FRG equation H87|l we see that and G and H as single frequency m = 1 terms will be fed by 0(?7A^) and 
O(^A^) respectively. Hence, given the rapid growth of 77 with scale, we expect these functions to be at least growing 
as rj with scale. We thus defined rescaled functions: 



Gi{u) = r)i^^ G(ue-C') 



Ad 

Hi{ui,U2,U3) = r]i 



(89) 



A4-2<i 2C/ ^ 



in terms of which one finds the flow equation for rj and G{u): 



9,^=(G'(0)-A"(0))77 

diG = (-2 + e - C)G + CuduG - 2A"G + (A(0) - A)G" - 3A'G" - G'(0)G - G"(0)A' 
+A'(2A"(0) + 2A") + 5iio(0, 0, u) + i(i?oio(w, 0, 0) - 2i/ooi(0, 0) - i/ioo(0, u, 0)), 



(90) 
(91) 



r 



Because of the above rescaling (|89|) no explicit rj appear 
in ^ 

Since G'(0) appears on the same footing as A"(0) ^ 
1/T/ in H9U|I it is natural to expect it to grow unbound- 
edly with scale in the same fashion. Indeed inspection 
of (|91|l reveals that G is fed by a term explicitly propor- 
tional to A"(0) which can consistently be balanced by 
the G'(0) term appearing in the first line of the same 
equation. Therefore we are led to expect that G itself, 
like A should exhibit a thermal boundary layer and the 



effect of temperature will be essential in understanding 
the structure properly. We will come back after taking a 
brief look at the m = 2 terms. 



2. terms proportional to square frequency m — 2 

The m = 2 terms, restricting to n < 2 -th order cumu- 
lant read: 



S. 



(1-2) 



iUrtDilrt-- 



iUrti *^irt2 [Urti Urt2^(Urti-Mrt2 ) + (Urtl-^^rt2 )^("rti^rt2 ) + ("rti-Urt2 )C'('"rti^rt2 )] (92) 



The renormalisation of D and of the functions A,B 
and C via a one loop calculation is performed in the 



Appendix. It turns out to be convenient to define: 

B,{u) = B{u) - ^C'{u) (93) 
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which simphfies the equations, the physics being ex- 
plained below. We define rescaled quantities as follows: 



FIG. 8: Shown is the G vertex (top image with no alphabetic 
label), and diagrammatic corrections to rj. Graphs (a-d) are 
contributions from tadpoles of the G vertex (note that (a) and 
(c) cancel by the same mechanism as dimensional reduction, 
and that (d) vanishes upon integration by parts on internal 
line). Graph (e) is the contribution from rjA. 



(k) (I) 



FIG. 9: Corrections to G. 



D = M^b 



A{u) 



Ad 



(94) 



Ad 



C(«) = ^C(^.e-C') 
Wi{ui,U2,U3) = rjf W{uie 



(95) 



U2e ^\u3e ^') 



and one finds: 



diD = (-2 - A"{Q))D - A{0) + C'{0) - ti^G'iO) - A"(0)) 

together with the FRG coupled flow equations for A{u),B{u) and C{u) which read: 

diA = (-d + Cua„-2A"(0)-4A")i-4A'i' + (A(0)- A)i" 

(2 G" (0) G' + 5 G"2 - 4G" A" (0) - 8G' A" + 2 A"^ + 4 G G" - 4A' G" 



(96) 



(97) 



diBi = {-d + Qudu - 3(A" + A"(0)))Bi + 3A"i?i(0) - 4A'i?; + (A(0) - K)]3'( + i(0)A" 
+rf{ - G'{Q)G' - G'2 - GG" + 2G'A"(0) + G"(0)A" + G'A" + A'G" - A"(0)A") 



(98) 



diC = ( - C - d + C,udu - A"(0) - 2A")G + A'A(O) - A'G'(O) - 3A'G' + (A(0) - A) G" + 2D A' (A"(0) + A") 
+rj^{ -2GG' (0) + 4A' G" (0) - 2 G G" + 2 A' G" + 2G A" (0) - 3 A' A" (0) + 2G A" - 2 A' A") (99) 

I 



From these equations we expect exponential growth of 
D, A, B, C at least as fast as rj^ due to the feeding terms. 
We expect from the considerations of Section lTTll that the 
growth is actually faster. Indeed considering the A flow 
equation at the origin gives 9/^(0) — — 6A"(0)^(0) keep- 



ing the largest terms of order 1 /T; and neglecting feeding 
terms. Note the similarity to the result of Section III. 
While we expect this qualitative behaviour, the precise 
nature of the growth is more subtle due to the fact that 
at non-zero temperature A{0) no longer satisfies a closed 
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equation. We will discuss this in more details below. 

Let us close this Section by noticing that all non-trivial 
terms in the r.h.s. of the above set of FRG equations (/?- 
functions) for G, A, B and C cancel when one chooses: 



G{u) = K'{u) 
A{u) ^ -rfK"{u) 

c = bK'iu) 
B, = b\K"{u) 



(100) 
(101) 
(102) 

(103) 



H{ui,U2,U3) = 3S'ioo(ui,U2,U3) (104) 
W{ui,U2, Us) = -3S'iio(ui, U2, Us) (105) 

and furthermore the flow of the kinetic coefficients sim- 
plify into: 



9/77 = 
diD^O 



(106) 
(107) 



This remarkable property, which serves as a useful check 
on the RG equations, can be understood in terms of a 
simple integrable model which is studied in the Appendix 
mi (which has very different physics from the one studied 
here) . 



C. dynamical thermal boundary layer 

1. dynamical action at non-zero temperature and FDT 

At r > two new effects must be taken into account 
not present at T = 0. First, in addition to the kinetic 
coefficients studied above one must take into account a 
variety of thermal noise terms. In the dynamical action 
this corresponds to terms with two or more u field having 
the same time argument. Second, new thermal contrac- 
tions are possible using the non-zero correlation function 
{uu) of the Gaussian theory. We first focus on the for- 
mer, and discuss the additional contractions at the end 
of this subsubsection. The action takes the general form: 



-|-oo n 



lUrti • • • *"rt„ 



fc=0 i=l 



(108) 



with = and the r*^ > are integers. The 3^]^ are 
translationnally invariant functions. Compared to the 
T — action it has additional powers of Ui and pos- 
sibly their time derivatives (such vertices are shown in 
Fig. I10|l . There is a temperature homogeneity degree 
s — X]fc Si '"f such that the term is ~ T^. The stan- 
dard thermal noise corresponds to Sp^Q ^ = —r/T, with 
= 2SkoSii. 



In order to maintain the Fluctuation Dissipation The- 
orem (FDT) there are relations between these vertices. 
A useful symmetry which constrains the allowed form of 
these terms is: 



Urt Ur.-t 



(109) 
(110) 



(we mean u_t — u'{—t)). For actions with no explicit 
time dependence, such as considered here, one can then 
later make a change of variables t —^ —t in integrals over 
times. We apply this to the bare action (jlll - I12|l . Con- 
sider first the infinitesimal variation of the interaction 
term: 

SSint ^ ^ I I iUruUrt2^{UrU " "rij) + O(A^) 



\r / iUrtA2R'iUrU " ""rtj + 0{\^), (111) 



using A(u) = —R"{u), a consequence of potentiality. 
This integrates to a boundary term which is a function 
only of the coordinates u and corresponds to the energy 
difference between the initial and final times(see below). 
Hence the interaction term is invariant for an arbitrary 
(r-dependent) A^. Unfortunately, this large symmetry is 
quadratically broken by Hll|l . First, the variation of the 
elastic term vanishes (up to boundary terms) only for 
spatially constant A,- = A. Thus the full action for = 
has a continuous global A-symmetry ()109f) . This can be 
used e.g. to put constraints on the terms appearing in the 
FRG equation order by order in 77 4i More importantly, 
the remaining terms in the action are only invariant un- 
der a discrete transformation, specifically (|109|l with: 



(112) 



Note that they are, however, exactly invariant (no bound- 
ary term) . 

Having established the invariance of the bare model 
under the symmetry (|112(l we know that it should be 
preserved under renormalisation. We must thus under- 
stand the consequences of this symmetry for a more gen- 
eral effective action. The relation to FDT is apparent 
since, performing the transformation in the path integral 
defining the response function one finds: 



Rt2-ti = (iut^ut^) 

= ((m_ti - u„t J-it-f^) = Rti-t2 + 



Cti~t2 

T 



(113) 
(114) 



i.e. the FDT relation for two point functions. The same 
is obtained from considering the action of the symmetry 
(|112|l on a general form (i.e. non-local in time) for the 
quadratic part of the effective action functional4^ 

We now discuss more precisely the conditions on the 
boundary terms which relate this symmetry to the FDT. 
This is simplest to see first in the context of the theory 
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before averaging over disorder. Let us define the Ito path 
integral: 

/•"(*/)="/ 

Z{uf,tf,u„U) = / DuDue-^"^' (115) 

J U{ti ) — Ui 



dufZ{uf,tf;Ui,ti) = 1 



(116) 



the (normahzed) conditional probability to find the sys- 
tem in state at time tj given that it is in state Ui 
at time ti. Here Sv is the MSR dynamical action in a 
given disorder realisation. By construction the Boltz- 
mann measure is the stationary distribution for this 
Z[uf,tf;Ui,ti) regarded as an evolution operator: 



dufZ{uf,tf;u,,t,)e-<^"<^^^^-''^^f^^/^ = 1(117) 



Note that if under the transformation above Sv Sv + 
SSv where SSv = ^Sv\u■^^ m/] is a function only of ui and 
u/ (boundary term) one has: 



Z{uf,tf;ui,ti) = Z{ui,tf;uf,ti)e 



-SSv 



(118) 



since t is changed to —t and thus boundary times ti and 
tf must be interchanged. Interchanging ut and w/ in 
(|117|l and using f lll8|) the normalization condition Hll()|) 
is found to hold only if: 

SSvKuf] - ^{H{ut^)-H{ut^)) (119) 

If on the other hand SSv depends also on the time deriva- 
tives at the boundary, then the FDT may or may not be 
satisfied depending on the precise nature of the boundary 
terms 4i ^ 

Upon averaging over disorder e^^~^^^^ — e^'^^^ , the 
shift SS obtained after transformation on the disorder 
averaged MSR action. It is readily seen that for the bare 
action, 6S is a sum of a single time integral cross correla- 
tion boundary term (and one response field) and a term 
with no time integral representing the second cumulant 
of the random portential V{u). More generally if one 
performs the transformation H1U9|) on the coarse grained 
model, one must obtain a SS which is a sum of bound- 
ary terms, each containing less response fields than time 
integrals. The n-th cumulant of the renormalized static 
disorder can then be retrieved from the corresponding 
boundary term with no time integral and 1/T" factor. 

We are now prepared to discuss one how can construct 
the effective action at finite temperature taking into ac- 
count the constraints from the FDT. It is useful to note 
that from the fundamental fields iii and u two linear com- 
binations transform simply under the symmetry (|112ll 



Y = 2Tiu - 



—u 



(120) 
(121) 



Terms in the effective action which are exactly invari- 
ant (i.e. whose variation do not produce any boundary 




O 



t < I 



FIG. 10: Compact notation for a generic vertex at T > 0. 
(a) 9 (b) „--f 



(c) 



(e) 



O O 



(d) Q 



- - - O 



-o- 



FIG. 11: One loop diagrams which correct the effective action 
at T > (in addition to the one for T = Q): the internal 
lines contain the full response function and the graphs are IP 
irreducible, (a -b) are the tadpoles using the full correlation 
function (the only possible ones as (d) should not be counted 
as it is a two loop diagram), (c) is the generic new one loop 
diagram at T > with two vertices (e) is one example of 
expanding the full correlation 



terms) must involve u only in the combination 2Tiu — ii. 
Examples will be constructed below. 

It is clear from these considerations that non-zero T 
terms can have time derivatives replaced by Tu. There- 
fore it is natural to group terms which formerly (at 
T — Q) were organized by the frequency power m by 
the more general index: 



M = 



(122) 



where Nu is the number of u fields appearing in the term 
and n the number of independent times, as before. Terms 
with a given Af , n can mix under the FDT transformation 

Let us start with M = 1 and n — 2. The only possible 
combination of the above invariants is, using symbolic 
notations as above: 



'T.1 



{2Tiui 



Ul 



8T 

-iuiiu2[{ui - U2) 



-J^iu^Giui ^ U2) + (1^2) 

T{iui - iu2)]G{ui - U2) 

(123) 
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recovering the zero temperature G term together with a 
non-zero T partner term. Note that at time ti we have 
used the invariant combination F^, taking care to remove 
the piece proportional to since there must be at least 
one response field at every time. This is not possible 
whenever there is an odd total number of fields ii and u 
at any particular time. Thus the iu2 above cannot be 
embedded in a term exactly invariant. Hence G{u) must 
be a total derivative and the above term gives a non- 
vanishing boundary variation under the transformation 
(|112|l . This term can also be understood before disorder 
averaging. It corresponds to replacing the naive zero T 
dynamical term ri{u, r)iuu corresponding to the damping 
in 1)82(1 by the invariant combination: 

/ (124) 

J rt 

Note that expanding out the factor in this term demon- 



strates that including u dependence in the damping coef- 
ficient has given rise to a white thermal noise which (for 
D{u,r) — g{u,r) — ) has w-dependent variance: 

(C(r, t)C{r',t')) = 2'n{u, r)TS{t ~ t')5{r - r') (125) 

The fact that G{u) is a total derivative then follows sim- 
ply from its interpretation as a cross cumulant of ri{u, r) 
and the conservative random pinning force /(m, r). Based 
on this reasoning it is clear that the function B{u) and 
G{u) must also be total derivatives. 

Note that the finite T partner of the G term is gener- 
ated in parallel to it from graphs of the form (e) in Fig. 1111 
One easily checks that it corrects the temperature term 
by SrjTuu where Sr] corresponds to the correction l|9U|) 
from G'(0) so as to maintain FDT relation. 

Let us now examine the terms M = 2 and n = 2. 
One can write the possible terms in a way which makes 
apparent the invariants: 



(2) _ 1 {ill - (2rmi - ui)^) {ill - i'2Tiu2 ~ U2?) 



T,2 - 2 4r 
(2rwi - u{)u\ ~ {2Tiui - iiif 
ST 

{2Tiui ~ iii)iil . ^ C"(ui — M2) {2TiiLi — ui)ui 



AT 

iu2Bi{ui - U2) + (1 



A{ui - U2) 
2) 



4T 



-IU2- 



in a way such that unwanted terms (with no u field associ- 
ated to a time) cancel explicitly, apart from the last terms 
where they combine to gives rise to a boundary term 
i9tj(ufC(ui — U2)). Accordingly, the B has been split- 
ted into B{u) = Bi{u) + G'{u)/2. The function Bi{u) 
must be a total derivative (see above) and its variation 
yields a boundary term. However, the invariance of the 
part cubic in the field in the Bi term is exact, which can 



4T 



iii2C{ui — U2) -|- (1 ^ 2) 



(126) 
(127) 
(128) 



r 



be traced to an exactly invariant term in the unaveraged 
dynamical action: 



{2TiUrt - Urt)ult ~ {^TiUrt - UrtY 

4T 



5(Mrt,r)(129) 



Expanding all terms above one can write explicitly: 



= ^tIo,2 + -^((*Wl)^W2W2 + (W2)^iwi'"l)^(ui - "2) ^ («Ul)^ («'«2)^^('«1 - U2) (130) 

3 

+ -T{{iiii)'^iu2Ui + {iii2fiuiii2)Bi{ui - U2) - T^((toi)''^W2 + iui{iu2f')Bi{ui ~ U2) (131) 
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Note that C does not give any bulk contribution at non- 
zero temperature. The FDT constraint only requires 
some boundary noise term for C. This is because C 
alone, with A = Bi = Q corresponds to a conservative 
dynamics 4i 

These considerations show that to the order studied 
the T > dynamical action is fully specified by the T = 



action and the FDT constraints. Thus we do not need 
to introduce at this stage any new operator associated to 
finite T. 

Having established that we are working with an appro- 
priate action (and hence have not neglected any pertinent 
coupling constants/functions), we turn briefly to the ef- 
fects of additional thermal contractions upon the FRG 
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equations. Up to this point, the only such contractions 
we have included are the "diffusion" terms (T/A" etc.) 
in the FRG equations for each coupling function. It ap- 
pears natural to neglect most effects of temperature since 
it is an irrelevant variable under the FRG, while clearly 
these diffusion terms are crucial, since they are necessary 
to stabilize the boundary layer. Within this treatment, 
the zero temperature rule of conservation of powers of 
frequency still holds. More generally, however, one can a 
priori perform thermal contractions that feed downward 
(i.e. reduce the number of time derivatives) in the fre- 
quency hierarchy, in particular by thermally contracting 
fields containing time derivatives. For some simple such 
contractions, a preliminary calculation shows that a can- 
cellation in fact occurs amongst the different "partners" 
required by the FDT, eliminating the unwanted mixing. 
We do not have, however, a general argument for such 
a mechanism of cancellation. Due to the complications 
of such a more general analysis, we will however proceed 
assuming this is generally true. We comment briefly fur- 
ther on this in the conclusion. 



D. dynamical boundary layer analysis: terms 
associated with averaged relaxation time 

Having established that to this order no new terms 
arise due to temperature we now attempt to study the 
structure of the thermal boundary layer in the operators 
studied so far. 

We consider the level m — 1 in some detail. We add 
the effect of temperature to lowest naive order which is 
to add to the right hand side of the term fG"{u), 
originating from the simple tadpole contraction of the G 
vertex. 

From examination of this equation we expect that G{u) 



supports a thermal boundary layer form for u ^ Ti/e 



(132) 



with g(0) = 0, g{x) an analytic function at a; = 0, odd 
and positive for a; > 0. It should match the fixed point 
form outside the boundary layer. For u ~ 0(1) and <C 
e we expect A(u),G(it) < A"(0) - G'(0) - e^xVf/. 
Thus in the outer region only three out of the several 
terms involving A and G are non-negligible. For now 
we neglect feeding from third cumulants functions S and 
H, we return to them below. The fixed point for G{u) 
outside the TBL is then trivially: 



G»=ra^-1 A'*(.) 

^ ' \ G'(o) ; ^ ' 



(133) 



Since at small argument A'*(m) — —ex it follows that 



g{x) g+a 



1 



2A"(0) 
G'(0) 



a constant, for large x. Using 



the TBL form to evaluate G"(0) yields: 

2 



5+00 = 1 



5'(0) 



(134) 



To analyze the boundary layer equation, we use the 
form H7U|) for A and similar forms for the third cumulant 
functions (|72|l for S and: 



H{ui,U2,U3) ^ (ex) h{ui,U2,U3) (135) 

fi 

The TBL equation for g{x) is then found to be: 



(136) 



= 2/".g + 3.97' + g'mf -g) + g"{i + /) + 2/'(/"(o) + /") (137) 

+siio(0, 0, u) + i(/ioio(", 0, 0) - 2/iooi(0, u, 0) - /iioo(0, u, 0)), (138) 
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where 5(711,^2,^3) = sj^ (ui,U2,M3), all these terms be- 
ing multiplied by e"^x/T/ while the terms originating from 
rescaling are proportional to e. This form will thus hold 
at scales such that T/ ^ e^. 

For given functions /, s, ft, this equation is an eigen- 
value problem for determining g'{0). This can be seen 
since for large x the linear problem has one exponen- 
tially growing solution, in addition to the one matching 
the outer solution which converges to a constant. Thus 
g'{0) must be tuned to select that solution. We illus- 
trate this behaviour in the approximation of neglecting 
all third cumulant functions. Then we recall that / sat- 



isfies: 

/" + /"(I + ./) = 1- (139) 

whose solution is f{x) — vT+~?— 1. Since, as discussed 
earlier G is a total derivative, it is possible to integrate 
the boundary layer equation once, and defining g = — /'+ 
7' one obtains: 

(1 + /)7" + 2/'7' - 5'(0)7 + (2/ - 1)(1 + g'm = 

with 7(0) = 7'(0) = 0. One interesting solution but unre- 
lated to the physics of interest here is 7 = 0, g'{0) = — 1, 
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i.e. g = — /'. It corresponds to an integrable set of 
models with a single exponential relaxation, which ex- 
actly obey the full FRG equations to one loop, and is 
studied in Appendix. A shooting procedure gives a solu- 
tion g{x) satisfying the proper boundary conditions with 
5'(0) = 2.646 ±0.001. 

Thus we find that the growth of rjl is determined by: 

9ir^=^(/"(0)+5'(0))r?7 (140) 
J-i 

yielding: 

r?7~exp(a(l)^) (141) 
Clearly a{l) = 3.646 is a non-trivial number. 



E. terms associated with second moment of 
relaxation time 

We now turn to the consideration of terms with m — 2. 
As emphasized in Section II the principal quantities of 
interest are the cumulants of the friction, the second one 
being embodied in A{u). The quantities Bi and C also 
appear at this order, complicating the analysis. Since 
these embody somewhat different physics we will focus 
initially on A{u) which fortunately satisfies equation H97|l 
which is independent of Bi and C . 

We add the effect of temperature to lowest naive or- 
der which is to add to the right hand side of H97I98I99|I 
the terms fA"{u), fB'{{u), fC"{u), respectively. These 
originate from the simple tadpole contractions. 

1. second moment of relaxation time: A{u) 

In the previous Section III we pointed out the rapid di- 
vergence of the moments of the friction (relaxation time) 
77,77'^' = ^(0), • • • driven by the low temperature diver- 
gence of A"(0). In doing so we neglected all functional 
dependence (such as A{u)). In the previous subsection we 
reconsidered the growth of the average friction rj^ which 
clearly does not itself has any functional dependence. In- 
stead the deviations of its growth from the prediction 
of Section III arise from a secondary mechanism of the 
feedback of G'(0) into rj. Physically it corresponds to the 
cross correlation G of the friction rj{u^ r) with the random 
force f(u,r) producing an increased growth ofrj. 

We would now like to reconsider the growth of the 
second moment A(0) = rj'^'^^ including functional depen- 
dence. In this case already the leading effect of enhance- 
ment due to the divergence of A"(0) is non-trivial. Thus 
we will focus on it here primarily ignoring secondary ef- 
fects of cross correlations between the random friction 
coefficient and the random force. In general these cross 
correlation effects enter the flow of A through G, H and 



W defined in H81|l . Terms involving H and W have al- 
ready been dropped in (|97|l for A. We will initially keep 
terms involving G in (|97|l but will drop them at a later 
stage of the analysis. It is not clear at this stage whether 
keeping these terms without simultaneously including the 
ones due to H and W would be consistent. 

We then note that A{u) satisfies a closed equation once 
G{u) is known. We first consider the nature of its solu- 
tion for u ^ 1 outside the TBL. Doing so one notes the 
presence of several large terms proportional to A"(0), 
G"(0). Balancing these large terms, we obtain the solu- 
tion outside the TBL: 

Ai{u)^tg^G*' ^-rfigl.K*" (142) 

where g^o was defined above. Note that, as was the case 
for G the convergence is very rapid due to the homoge- 
neous part dA = —2G'{Q)A. The important feature of 
this result is that A{u) is of order rjf outside the TBL. 

We are going to search for a TBL solution for A which 
grows faster: 

A{u)^t'^hC-^) (143) 

with A > 2. In order to match the above solution outside 
the TBL one should have h{oo) = 0. The TBL equation 
for h then reads: 

(-A(/"(0)+g'(0)) + 2/"(0)+4/")/i 

+ {1 + f)h" + 4f'h' = (144) 

Due to the presumed faster growth of A than rjf (A > 2) 
the feeding terms in (|97|l are negligible and have been 
dropped. As discussed above, since our principal inter- 
est is to compare the growth of the second relaxation 
moment parameterized by A{0) relative to the growth 
of the mean rj, to be consistent we drop the analogous 
renormalization oirj by G"(0), i.e. set ^'(O) = in (|144|l . 
Numerical solution then yields: 

A = 2.64.. (145) 

The analysis is thus consistent since we find A > 2. To 
the order considered we therefore have: 

V^i-Vf''-»VI (146) 

This gives, in the present approximation q;(2)/q;(1) — 
2.64. As seen in the previous subsection we expect both 
a{2) and q;(1) to be both increased by inclusion of the 
effect of cross-correlations of friction and random force. 

2. growth of other 0{uj^) kinetic coefficients: D, B, C 

Due to the feedback of A{0) into D we expect D to 
grow at least as fast as rjf. In the simplest scenario, 
indeed, all quantities would scale the same in the same 
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manner. However, we see no general reason why this need 
be the case. Indeed, examination of Bi and C using the 
same truncation scheme as for A, shows that they grow 
faster. We sketch this analysis here. Consider first C{u), 
which also feeds into the inertial mass D. We assume C ~ 
rj'^, with fi > \ > 2. With such growth of C, the feedback 
of C"(0) into D will overwhelm all other feeding terms, 
and we expect D = rj'^D, with D scale independent. It 
is then natural to define C{u) = —DC{u). H96|) becomes 



Thus to leading order in l/T), one needs 



C (0) = m(G'(0) - A"(0)) 



(148) 



diD = (-2 + ^i{G'{{)) - A"(0)) - C {0))D. (147) Using the above forms, C satisfies 



dC = (2-d + C9«- A"(0) + Ai(A"(0)-G"(0))-2A")C + fiC 
-A'C'(O) - 3A'C' + (A(0) - A) c" - 2 A' A"(0) - 2A' A" 

I 



(149) 



to leading order, i.e. dropping terms ~ 77^"^, and 
neglecting feedback from higher cumulants as before. As 
for G and A, the outer solution for u ~ 0(1) is readily 
found equating the large terms ^ A"(0) + C'{Q) ^ 1/fi: 



C 



C'jO) + 2A"(0) 
(/i-l)A"(0)-/.G'(0) 



A'(u) u-O(l). 



(150) 

As before, for small u we make a TBL ansatz, 

Ciu) = exc{exu/fi), (151) 
which yields an equation very similar to (|138(l for g{x): 



(1 + /)c" + 3/'c' + (/"(O) - ^iif"{0) + g'm + 2f")c - /'(2/"(0) - c'(0) + 2/") = 0. 



(152) 



r 



We require, to match (|150|1 . that c goes at a constant at 
large argument, and c(0) = since c is an odd function. 
Furthermore, from ifTtS)) . we have c'(0) = -/x(/"(0) + 
(/'(O)). This formulates an eigenvalue problem for ^. As 
above, to solve, we use the (approximate) form for f{x) 
in H139|l and, for consistency as before set g'{Q) = Q. A 
shooting procedure gives ^ — 3.377, indeed greater than 
A as required for consistency. In summary we find the 
growth of the kinetic coefhcients: 



C{u) - b 



, 7^3.377 



(153) 



Finally, we discuss the growth of Bi . Since it is fed by 
A{0), it must grow at least as fast as rj^, so all other feed- 
ing terms on the last line of H98|l are certainly negligible. 
Remarkably, even in the presence of the thermal TiB'^ 
term an asymptotically (for large I) exact solution can be 
found. In particular, one finds that the homogeneous (in 
Bi) part of the Bi equation has an exact eigenfunction 
which is just Bi{u, I) = Bi{l), a constant independent of 
u. This turns out to be the most unstable eigenfunction, 
with eigenvalue —d— 3A"(0). The exponential growth 
of this unstable eigenfunction is faster than that of A{Q), 



and hence dominates the flow at large scales. Hence, 
writing this relative in terms of rji (neglecting the g'{0) 
renormalization of 77 as before), one finds 



Bi{u) 



(154) 



CONCLUSION 



We have through a series of successively better approx- 
imations arrived at a description of the growth of the mo- 
ments of relaxation times (friction coefficient) encoded as 
eigenvalues of functional FRG equations. This final stage 
of analysis was carried out only for the mean and variance 
- the extension to higher moments is a formidable techni- 
cal challenge. Nevertheless, already at this level we have 
observed how these functional eigenvalue problems pro- 
vide a mechanism for describing a broad but non-trivial 
(i.e. not log-normal) distribution of time scales. This 
is at variance with numerous other existing examples of 
systems exhibiting simpler log-normal tails which can be 
obtained from simpler non-linear sigma model diagram- 
matic calculations, such as in disordered conductors jM- A 
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similar log-normal tail was indeed obtained in Section III 
from an approximate truncation of the FRG equation. A 
rather strong physical difference from the aforementioned 
quantum diffusion problem is the rapid exponential scale 
dependence of the relaxation times for 9 > 0, very dif- 
ferent from the logarithmic dependence of two dimen- 
sional weak localization corrections. It is an open ques- 
tion whether some less trivial distribution might arise at 
the metal-insulator transition in c? > 2 and whether sim- 
ilar functional renormalization ideas might be useful in 
this context. 

Many outstanding issues and extensions remain. Of 
these, the most fundamental are germane to both the 
dynamics and the statics. In particular the very ba- 
sic problem of perturbative control of the theory (most 
interestingly in the e-expansion) remains unsolved. This 
question, and the associated matching problem of relat- 
ing e.g. random force quantities like the fk in l|72() de- 
fined deep within the boundary layer at w = to the zero 
temperature ones occurring far outside for u ~ 0(1) are 
better addressed in the simpler context of the statics. We 
will refrain from commenting further upon them here. 

Of the problems specific to the dynamics, perhaps most 
important is a systematic treatment of all thermal terms 
in the FRG. We have begun this program by classify- 
ing all operators associated with "thermal noises" in the 
effective action (Sec. liV C 1|) consistent with the FDT. 
However, up to this point we have included the effects 
of non-zero temperature only through the leading "diffu- 
sion" terms (T;A" etc.) in the FRG equations for each 
coupling function. As discussed in Sec. lIVCTTl while this 
assumption is natural, we do not at this stage have a 
general justification for it. The importance of additional 
thermal contractions thus remains an important issue for 
further investigation. 

Once these basic remaining issues in the FRG formula- 
tion are resolved, the present methods offer the opportu- 
nity to explore numerous physical problems. Obviously, 
equilibrium response and correlation functions are of con- 
siderable interest. Perhaps the approximate techniques 
of Section III (and Appendix IE|) may have an extension 
to the full functional description. It will also be valuable 
to reinvestigate the response to a uniform applied force 
in the creep regimCj^?. in light of the full dynamical struc- 
ture of the thermal boundary layer exposed here. Appli- 
cations of these ideas to non-equilibrium response and 
aging is also tantalizing. Similar approaches should be 
applicable to quantum problems in the Keldysh formal- 
ism. These and other applications of the present formu- 
lation certainly provide a broad scope for future progress 
in understanding glassy dynamics. 
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APPENDIX A: SINGLE TIME SCALE 
CALCULATIONS: EQUILIBRIUM 

1. analytical results for the equilibrium response 
function 

It is interesting to observe how the two putative scal- 
ing regimes described in Sec. ^ arise in a detailed calcu- 
lation of the mean response function. To do so, we de- 
velop an FRG scheme to calculate directly the response 
function at arbitrary uj, k within the scaling regimes de- 
scribed above. It is necessary to follow the flow of the full 
wavevector and frequency dependence of the "kinetic" 
part of the MSR action. We therefore generalize the form 
in ini) to: 

Sq[u,u] = / iUrt{Rf^)rt,r't'Ur't' 
J r.r' ,tt' 

-r]T I {iurt){iurt), (Al) 

Jr,t 

where, in a slight abuse of notation, we have denoted the 
quadratic MSR kernel by R^^. Using I|A1|I . we extend 
the FRG analysis leading to IjlSI) to derive an RG equa- 
tion for the response function. As before, the strategy 
is to integrate out spatial Fourier modes A > k > Ae~\ 
but now keeping the explicit time dependence. At this 
stage, we will not assume time-translational invariance, 
though we will specialize to this at a later stage. The 
FRG equation for R^^ is 

diR-j{t,t') = -TiAf{RA„i{t,t')-St,t'J dt"Rf,^,i{t,t")) 

(A2) 

where ti is an initial time at which the system is prepared 
in some as yet unspecified state (or distribution of states). 
(|A2|) is obtained formally by computing the correction to 
the (inverse) response function upon integrating out the 
modes in the shell, and using definitions IjlSI) and 121|l . 
We perform this integration perturbatively in A (to first 
order), which gives the lowest order term in e. 

At the end we want the true response function Rq{t, t'). 
It will be obtained by integrating the flow from I = with 
the initial condition: 

R,,iMt,t')^e-'^'^'-''^eit-t') (A3) 

setting — 1 for convenience, up to the scale I* such 
that Ae^' — q. 

Rq{t,t') ^ RqJ = U.{A/q){t,t') (A4) 
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This amounts to neglect contributions coming from the 
modes k < q, a.s is usuaUy done in the RG. These are 
examined below. 

Although the initial condition in (|A3|) is time- 
translationally invariant (TTI), the solution of the RG 
equation does not in general remain so, due to the pres- 
ence of the initial time t^. This leads to the aging proper- 
ties to be discussed in the Appendix IbI The TTI regime 
is recovered in the limit ti ^ oo (for large but fixed fi- 
nite size system) where one can set Rk{t,t') = Rk{t — t'). 
Then (jA2|) can be Fourier transformed in t — t', Rkiioj) = 
/t>o Laplace transformed with s = ius) to 

obtain 



Rf. ^{ioj) ~ iuj + k"^ + 



1 



1 



(A5) 
:) 



(A6) 

where we have defined a "self-energy" 'E,k{iuj) with ini- 
tial condition 'Sk=i{ii^) ~ 0. To obtain (|A6p one writes 

SA;(icj) = /o"^'^^'^' diR'j^\{iuj), uses the Fourier tranform 



of (|A2|) and differentiate w.r.t. k (we set from now on 
A = I). One can check that consistently Efc(O) = 0, as 
requested by the statistical tilt symmetry, which we use 
from now on. Apart special cases^ (|A6p does not ad- 
mit analytical solution and we now analyze the various 
regimes of interest. 

From ljA6p one first finds the small lo behaviour of 
RZ^iiuj) as: 



i?,^, ^{iuj) = k'^ + iujrjk + 0{uj'^) 



(A7) 



T.k{iLu) ^ k'^g{y), 



(A8) 



in (|A6|I gives the closed differential equation yg' — g/{l + 
g), which has the implicit solution (taking into account 
the behavior of S for iuj ^ from ljA6|l ') 



(A9) 



(|A8IA9|) correspond to the Y scaling limit of Sec. ITU 

(jASp is valid for finite y. As y ^ cxd, we enter the 
logarithmic {X scaling limit) scaling regime, in which 
the scaling variable x = fc(ln(l/itt')//3)^/^ is fixed and 

cxD in this limit, the first term 
U can be neglected, leading 



iuj, fc^ ^ 1. Since g{y) 



on the right hand side of 
to the ansatz 



(AlO) 



with (2 — 0)f + xf = — 1 from (jASp which determines 
the form of the scaling function f{x) as: 



1 



2-e 



1 



(All) 



and the constant Xc = {1/6)^/^ is determined by match- 
ing to (jA8|) . This regimes exists only for: 



k < 2;c(ln(l/it^)//3)-i/^ 



(A12) 



and thus corresponds to the limit of small wavevectors at 
fixed w, or to relaxation times t <^ Tk (r ~ e'^*^^/'^) for 
X < Xc)- When x ^ x~ one crosses over to the y regime. 

We can check that these results merge smoothly with 
the result directly obtained at the upper critical dimen- 
sion c? = 4. There the equation for S]fe(iw) becomes: 



1 



(ln(I//e))^ uv + k"^ + hk[iuj) 



1) 
(AI3) 



which yields the same two scaling regimes, the first one 
with Tfe ~ /c~^ exp(/3/(2fc^(ln(l/fc))^)) and the same scal- 
ing function g{y) (|A9p and the second one reading: 



s(fc) = 



13 



ln(l/fc) 



(-1 



ln(l/fc) 



ln(ln(l/itj)//3) 



We now turn to the calculation of the response function 
in the time domain. Consider first the regime y — t/rk 
fixed and t —^ oo, k ^ 1. Inverse Laplace- Fourier trans- 
forming l|A5|l and using H24|l gives the scaling form 



where rjk satisfies dkTjk — —Pk^^^^rjk which yields rjk — 

k'^Tk, i.e one recovers as expected the single caracteristic -^^rith 

time scale given by if^ . 

To analyze further the higher order terms in iuj from 
(|A6I) . we first consider the scaling regime iw, fc^ <C 1 with 
y = iujTk fixed (which implies iu <^ k"^). Making the 
scaling ansatz 



Rk{t) 



1 



A;2 



Tk 



200 + 7 



— 200 + 7 



i + 5(y)' 



(A15) 



(A16) 



While we have not performed a complete analysis of the 
integral in (|A16|I . the large time behavior can be ex- 
tracted .-^ For 3^ ^ 1, the integral is dominated by the 
vicinity of the branch point on the real negative axis at 
y — — 1/e, leading to 



Rk{t) 



1 



Tk 





t 








eTk_ 



t-»Tk. (A17) 



In the logarithmic scaling regime, we cannot simply 
invert the Fourier space result in (|Af 0|l . as it does not 
extend over the entire frequency interval. Instead, we re- 
turn to the defining RG equation for R'i^^{t), ljA2|l . By in- 
verting this formal equation, and again integrating down 
to scale fc, we obtain an equation for Rk{t) directly: 



dkRkit) = -2kRk *t Rk 
-~pk^-\Rk H Rk H Rk - Rk *t Rk 



(AI8) 
Rkit')) 



t'>0 
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where *r denotes a convolution. Note that, aside from 
the momentum-dependence of the prefactors and the ab- 
sence of derivative terms, (|A18(I bears a formal similarity 
to the mode-coupling equations of mean field models. In 
order to match the scaling expected from the above log- 
arithmic frequency regime, we make the ansatz 

"''"'- .(i„.)^-V.g./. ^i'^4''°'°" 

with F[0] a constant. Inserting this in IjAlSp . it is 
permissible to drop the first two terms in the loga- 
rthmic scaling regime, and moreover to approximate 
dt'R{t - t')R{t') w 2R{t) jl dt'R{t'). This yields: 

eXF'[X]^2F[X] dzz-^'^-^/'^^F[z]. (A20) 
The solution is 



-f-oo 



dh 



in{x^-T^ F[x]) ^{2-eY + A6eh ye) ^'^('^/'^ )' 

(A21) 

where X* — 1/6 in the boundary of the regime, at which 
F{X) diverges, signaling the onset of a regime of faster 
relaxation onto the regime y. 

Having computed the response function Rk{t) in the 
equilibrium TTI regime, we also obtain the time depen- 
dence of the connected correlation defined as Ck{t) = 
(Mfe(t)u_fc(0)> - (ufe(t))(u_fc(0)), with Ck{t ^ +(^) = 
and Ck{t — 0) the equilibrium connected correlation. 
Indeed they are simply related through the fiuctuation 
dissipation relation dtCk{t) = —TRk{t) or, in frequency 
space, 



Ck{iuj) = ^{Rk{iuj) - Rk{-iuj)) 

lU) 



discussion 



(A22) 



We now pause and comment on the results of the FRG 
calculations we have just obtained. Let us first men- 
tion the nice features before stressing the unsatisfactory 
points below. 

First we note that the Wilson scheme used directly 
on the response function within the single time scale as- 
sumption indeed yields, as we anticipated from general 
arguments in Sec. ^ two distinct scaling regimes, the 
y = t/Tk regime and the X = k^lnt regime, with scal- 
ing forms in HA8IA10IA15IA19|) . The scaling functions 
themselves were found to be non-trivial, with interesting 
analytical structure. While the existence of the y — t/rk 
regime seems to be a straightforward consequence of the 
assumption of a single time scale , the emergence of the 
regime within this hypothesis is less obvious. 
Within the Wilson scheme, it seems to result from the 
system keeping a memory of a whole spectrum of smaller 
relaxation times r ~ (■7fe)^/^'=, x < Xc, generated during 



the coarse graining procedure and naturally appears here 
(while one would naively expect the largest one only, to 
play a role). It does have the form of activated dynamics 
since the scaling variable is truly X — {T/T*)k^ \nt and 
thus corresponds to crossing barriers of heigth ~ k~^ . 
That such an activated regime should exist is physically 
rather natural. Indeed we expect from simple droplet 
arguments that the equilibrium dynamics of mode k at 
large time difference (in general t — t' , denoted here t) 
is dominated by the rare active configurations with (at 
least) two quasi-degenerate low free energy states at scale 
L ~ 1/fc of the system. The probability to find two nearly 
degenerate minima (on the scale of the thermal energy T) 
at scale 1/fc is Tfc^. One expects these two minima to 
be separated by a barrier Ut also scaling like Ub ^ k~^ 
and thus when Tint > Ub equilibrium thermally acti- 
vated motion back and forth over this barries^ becomes 
active and gives rise to time-dependent correlations on 
the scale Int ~ fc~^. Our analytical result thus exhibits 
the correct scaling behaviour and it is thus encouraging 
that such barrier crossing behaviour and scaling comes 
out of the present RG calculation. 

Upon a closer look to our results in the X regime, ev- 
erything works as if there is an effective distribution of 
smaller barriers Ub = x'k^^ with a distribution of relax- 
ation times T = e^'^'' for < x' < Xc- The total weight 
of this distribution being only ~ Tfc^ it can be written 
as k^ / (34>{x'). 4>{x') diverges at x' = Xc- This is easily 
seen, e.g. on the form for the correlations. Indeed, using 
the above FDT relation one obtains: 



C'kiti) - Ck{t2) 



fc2 (i 



k" \nt2/0 



du 



F{u) (A23) 



'kO\nti/P U 

for the correlations in the logarithmic regime. In this 
expression, the T/fc^ equilibrium correlation is usually 
explained as T/P = (l/fc'^+2?)(rfc^), i.e the product of 
the size of a positional fluctuation between two degen- 
erate states at scale fc and the probability of this active 
configuration to occur. Thus we see that there is here 
an additional reduction by an extra factor Tk^ , the to- 
tal weight of barriers much smaller than Tk (note that 
the above correlation variations within regime X are sub- 
dominant compared to the ones in regime Y, which really 
account for all but a small fraction of the total variation) . 
Similarly one sees that the response corresponding to a 
barrier x'k~~^ can be written as: 

^ ^ ^ ^ \e{x-(A)2i) 



fc2 I 



fc2 1 ^ gPk-Oix'-x) 



the weight jk^(f>{x')dx' yields exactly our result Rk{iuj) 



as fc — > with X ~ k^ \n{l / ioj) / (3 . Thus averaging with 
)dx' yie 

in regime X if one chooses: 

= ^."^2-0 1 (A25) 



2- 



1 



Another puzzling feature of the above results is the 
nonmonotonicity of the above scaling functions. As dis- 
cussed below this is directly related to the assumption of 
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a single time scale, and has prompted us to reconsider 
the whole calculation (at a high price of technical dif- 
ficulty) in Section ITTll We see from l)A6|l that T,k{icu) 
is a decreasing function of k always and that Rk(ii-tj) 
at fixed iuj is an increasing function of k for k small 
enough. Correspondingly, the real-time solution Rk{t) 
is an increasing function of k at fixed t throughout the 
logarithmic regime and also in the short-time portion of 
the t ^ Tk regime. Similarly, ljA23|l implies that the 
correlations are also increasing functions of k at fixed t. 
While this behavior is unexpected, we are presently un- 
sure whether it is in fact unphysical. What is clear is 
that it is a consequence of the single time scale assump- 
tion. Indeed, the rather simple and apparently physical 
expression Rk{t) — e~*/'^'=/(fc^rfc) is also increasing with 
k for small t/r^. It appears that one can argue fairly 
generally that, provided there exists a long-time regime 
with a well-deGned r^, the response must be increasing 
with k for t/rk 1- 

As discussed at length in the text, one does expect 
that the single time scale description is unsufficient and 
one should instead consider a distribution of time scales. 
Let us examine the question of monotonicity when Rk (t) 
is simply a superposition of elementary relaxation pro- 
cesses. Discarding the subdominant prefactor and 
writing Tfc = e^^ where L = 1/fc is the scale, one can 
consider the average: 



Rkit) = J dUPLp] 



(A26) 



It is dominated by the saddle point U*{L,t) solution of: 



(A27) 



and the condition for Rk{t) to be an increasing function 
of L is 

dLPL[U*]>Q (A28) 

For a gaussian \nPL[U] = -{U - L^)V(2sL") and 
t ~ (the worse point) one finds U* — — sL" and 
\nPL[U*] -i* + fL". Thus one needs a > 9. Note 
however that this supposes the gaussian to hold down to 
U* < 0, which may not be the case in general. On the 
other hand the only real condition concerns the mono- 
tonicity of the scaling function itself. Thus, one way to 
reduce the effect of nonmonotonicity is to increase the 
width of the distribution of time scales. 

To close this discussion, it is useful to contrast the 
present situation of an elastic system with fast growing 
barriers with what happens in the marginal case 9 — 0. 
This is realized for a periodic model in d = 2, e.g. for the 
line of fixed points of the Cardy Ostlund model. There 
of course one expect a single scaling regime compatible 
with simple matching arguments. Setting 9 = in (|A6|) 
one finds the exact solution: 



fc2 = (l 



Sfe(iw) 



)-2//5(i + jc^ ^) (A29) 



2 + (3 



2 + P 



(A30) 



obtained writing -d{k'^)/dE = (2//3)(l -I- k'^iiuj + E)-i) 
For w <C 1 this yields the scaling form 



T,kiiuj) ^ k'^giy = iujk ^), z = 2 + /3 



y = .9(1 + 



2 + f3 



9) 



/3/2 



(A31) 



where z is the equilibrium dynamical exponent (note that 
for (3 -foo one recovers (|A9|) . The self energy nicely 
interpolates between Sfe(ia') ~ iujk~^ at small iuj <C k^ 
(as also obtained from considering the flow of the uniform 
77; ~ e'^') and Sfe(itj) - (f )^/^(iw)2/^ at large iuj > k\ 
^From there one obtains the response function Rk{t) — 
k^Q{tk~^)^ which is found to decay as in IIA17|l with Tk ^ 
k~^ and a characteristic time (1 -I- 2/j3)^/'^Tk (instead of 
erfe obtained for [3 — > +oo, and which behaves as: 



Rk{t)^r 



(A32) 



in the limit 1 < i < fc" 
equation: 



The function Q obeys the 



+ (2 + ~p)yg' = (-2 + 0)g *yg - pg*yg*yg 

(A33) 

Note that such as scaling function g of y = tk^ leads to 
a trivial scaling regime in X = lni/ln(l/A:) reduced to a 
delta function at X = —z. Note finally that even in this 
case, the scaling regime Rk{t) is again nonmonotonous: 
it vanishes at /c = 0, increases up to k* , with t{k*Y = zf3 
and decreases beyond. 



APPENDIX B: SINGLE TIME SCALE 
CALCULATIONS: NONEQUILIBRIUM AND 
AGING 

1. response function and various regimes 

The RG recursion relation for the response function 
derived above was not restricted to equilibrium, and 
it is thus interesting to write down the corresponding 
equations (for response and correlations) in the full non- 
equilibrium regime. Within the intrinsic limitations of 
the single time scale approach, this allows in principle to 
acces the aging properties of the system. 

To obtain closed equations for the two time response 
function Rk{t,t') one again iterates ljA2ll from the same 
(TTI) initial condition Rk,i=o{t,t') = 9{t - t')e-'='(*-*') 
up to ^ = In (1/(7), keeping finite and making no TTI 
assumption. Thus the response satisfies the differential 
equation: 



R-\t,t') ^ Su'idf + k'') + i:k{t,t') 



(Bl) 



dki:k{t,t') ^ f3k^-''{Rk{t,t') - Sw J dt"Rk{t,t")) 
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where matrix multiplication and inversion is with respect 
to (t,t'). It is more convenient to avoid the two time self 
energy and write a closed equation for Rk{t,t') as: 

dkRk{t,t')^-2k{Rk-Rk)it,t') (B2) 
-^k^-\iRk-Rk-Rk)it,t') 

- f dhRk{t,h)Rkih,t') f dt"Rk{h,t"))iB3) 
Jt' Jti 

with initial condition Rk=i{t,t') = e{t - i')e"'=' 
The full analysis of this equation is quite complicated 
and we have not attempted it. We will give only a few 
features, at a naive level, which remain to be confirmed 
by a more detailed analysis left for the future. 

The function Rk{t,t') depends on three variables but 
in the limit k ^ 1, t' — ti ^ 1, t — <' ^ 1 we expect 
that it takes scaling forms depending only on two vari- 
ables. What these variables really are depends on the 
time regime, and one can identify several possible time 
regimes and subregimes. They can be classified as fol- 
lows, where we indicate the form expected for the re- 
sponse function i?fe(t, i'), by order of increasing time and 
time differences: 



{I) <1, 

Inrfc 

Ht-t') g(fc'^ln(t^tO,Ml_p) 
^ ' Int' (t - t'Mn'' (t - 1') 



{lb) t~t' ^t' 



{t-t')lrP{t-t') 
h{k^ In(i-t'), V") 



Int' 



Int ln< /(fc^'lni, 

{Ic) 7-7 > 1, -j < 1 • 

In t' In Tk 

(Id) t^Tk-. 



t' In" t' 



(B4) 



{II) t'^Tk, 



Ht-f)^^, F{knn{t-t'),fj 



(Ila, — . . . j- 

liiTfc {t-t')ln''{t~t') 



{III) 



luTfe 



fc^Tfe 



>1 equilibrium Rk{t~t') 



(B5) 



{Ilia) Mi^<l F{knn{t^t')) 



In Tfc 



{t^t')\u^-^'\t~t') 



Q{—) 

{I lib) t^t'^T, 



Regime (III) is the equilibrium TTI regime, where the 
only dependence isint — t'. There are two scaling forms 
possible corresponding to the two subregimes X (Ilia) 



and Y (Illb) studied in Appendix 1X1 Fully equilibrated 
regime (III) is expected here for very large times t > 
t' S> Tfe, and is somehow at variance with mean field 
models (where one always expect aging, e.g. for t ~ 
t', even for very large t'). In regimes (I) and (II) the 
mode k at t' has not yet equilibrated, and the scaling 
functions ar now also function of k^ In t' (regime I) or 
t'/rfc (regime II), in addition of being functions oit—t'. In 
both regimes (I) and (II) if t— t' is small one expects some 
kind of equilibrium regime. Indeed for t — t' ~ 0(1) we 
expect that there will be a fully TTI equilibrium regime, 
but it is also expected to be nonuniversal. A universal, 
quasi- equilibrium regime is expected however for t — t' ^ 
t'" < t' , u < 1 (regimes (I a) and (II a)). As the time 
difference increases it should crossover at < — i' ~ i' to 
an intermediate aging regime (regimes (I b) and (II b)). 
Regime I is most complex as there one expects two later 
regimes as i — t' ~ t ~ t'", > 1 crossing over to yet 
another scaling regime when t reaches . It is interesting 
to note that either in Sinai model^ or even more clearly 
in the ID random field Ising modeli such regimes are 
also expected, some have been studied demonstrated and 
studied in details (there is also an equilibration time scale 
analogous to Tk). 

In the determination of all the above regimes the quan- 
tity: 



Mfc(^) 



dt"Rk{t,t") 



(B6) 



which appears in ljB3|) plays an important role. It is a 
function of t alone. It satisfies the equation: 



dkMt)= / dhRk{t,h)[~2kfik{ti) 



dt'Rk{h,t')nk{t'))] 



(B7) 



Pk'-''{Mti)' 



Its value can be understood by using the STS covariance 
under Urt — > Urt + Vr, where Vr is an arbitrary func- 
tion. In a general non-equilibrium situation, the STS 
gives constraints relating different initial conditions at 
t = ti. It can be written as hiZ[hkt, hkt,Uk.t=o = 0] = 
Z[hkt + k'^Vk, hkt, Uk,t=o = Vk] - hktVk, where the ini- 
tial condition is explicitly indicated. It thus immediately 
yields: 



S < Uq{t) > 
5Uq{ti) 



(B8) 



When ti ^ Tk we expect that the influence of the ini- 
tial condition on mode k has been washed out, and 
we find the equilibrium constraint limt^+oo M*:(i) = 
/o °° Rk{T)dT — 1/k^ (which, combined with FDT gives 
Cfe(O) — Cfc(oo) = T/k^). Thus, we expect fik{t) to take 
the form: 

fik{t) = k-^+''m{k'^lnt) < 1 (B9) 



lJ-k{t) 



luTfc 

k-^M{t/Tk) t^Tk 



(BIO) 
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with A^(oo) = 1 and a reduction k"' in the short time 
regime compared to asymptotic one, with an interpre- 
tation in terms of the susceptibihty to initial condition 
being almost 1, presumably with some rare (droplets ?) 
configurations exhibiting decorrelation. 

To discuss the specific choice of the scaling functions 
and prefactors we proceed as follows. Let us consider 
regime (I). We have found that with the forms of the 
prefactors in subregimes (I c) and (I a) indicated above 
we could obtain from HB7|I non-trivial equations for the 
scaling functions. The regime (I b) is then necessary 
to match (I c) and (I a). Next, with the forms con- 
jectured for (Ia,b,c) the terms in (|B7|1 scale respectively 
as 1, A:2-e(i-a0^ fc4-e-2e(i-a.)^ ^2-e-e(i-a0+^ (^here 

the first term is the derivative) with at — 7, (5, a re- 
spectively in each subregime. Thus either a < 9/2 and 
6{1 — a) = 2 — 9 + (J and only the last term counts 
or (J > 9/2 and 9{1 - a) = 2 - 9/2 (in equihbrium 
regime III one had a — leading to a = 2 — 9/2). 
Here, we see that if /ifc(0 is determined by an inte- 
gration over regimes (Ia,b,c) as is natural, it implies 
CT = 2 — maxQ,;=^^5_Q, 9{1 — a^), and one sees that a = 9/2 
and: 

7-<5 = a = ^-^ (Bll) 

is the only solution. Note that this contradicts the naive 
expectation that the form in the quasi equilibrium regime 
(I a) would scale as the equilibrium form Ilia (they differ 
by a power of a ln(i — t')): to get 7 = 2 — 2/9 would 
require cr = 0, and some argument that the value of fJ,k{t) 
is controlled by i — t" in the short time non-universal 
regime. 

Accepting the above scenario as reasonable we find the 
equation for the scaling function f{x,u) as: 

9xda:f{x,u) = /3a:2(i-") x ( (B12) 

dui U r' du2 , W2 

/ —JKX,ui)f{xui,—) / —f{xui,—) 

dui du2 U2 u 

~ — —j{x,Ui)f{xUi,—)f(xU2,—)) 
Ju Ui U2 

The bound in the integral really comes from the ration 
\n.ti/ \nt assumed to be very small. The first term in the 
right hand side of (|B3I) gives a subdominant contribution. 
Similar equations hold for the other regimes. We have 
not attempted to analyze further these equations at this 
stage. This would be necessary to fully confirm the self 
consistency of the scenario proposed here. 

2. correlation function 



full local quadratic term in the running effective action: 

- \ I {iUrt){lUrt')Ul{t,t') (B13) 

^ Jr,t>ti,t'>ti 

One can also decompose it as Ui{t,t') = Vi{t,t') + A;(0) 
by extracting the persistent part (disorder), requiring 
that \imt.t',t-t' ^+00 Vi{t,t') = 0. To lowest order 0(A) 
Ui is corrected and flows as follows: 

-CAe~',l{t,t')) (B14) 

where we assumed a flat initial condition Urt=a — (oth- 
erwise it should be added) and to this order Ui remains 
local. The persistent part of (|B14|) yields 9;A;(0) — 
-TTiAf in agreement with l(TK|l . (fTTjl (using that the per- 
sistent part of the parenthesis in (|B14|) is the equilibrium 
connected correlation C'^g-i ; = 7]A^^e^'). Substracting 
it yields the flow of V/. One closes the equations deter- 
mining Ci, Ui using Cfe,i(t,i') = [Rk,i ■ Ui ■ RkA{t,t'). 
Equivalently one can separate the effect of the ran- 
dom force part of the disorder in the correlation and 
write Ck,i{t,t') = Ck,i{t,t') + Ai{0)fXkAt)(^kAt') Mh 
fJ-kd{t) = Jq dt'Rk,i{t,t'j) write two closed equations for 
Ci and Vi using CkA{t,t') = [Rk,i ■ Vi ■ RkA{t,t'). 

Proceeding as above, to determine the correlation one 
defines Uk = Ui=\n{A/k) (similarly for Vk) and obtains: 

dkUkit,t') - pk^-'{^Ck{t,t) + ^Ck{t\t') - Ck.{t,t')) 

(B15) 

with Ck.i=in{A/k){t: t') = Ck{t, t'), which should be solved 
along with: 

Ck{t,t') = J dtij dt2Rk{tM)Uk{tiM)Rk{t' M) 

(B16) 

with initial conditions at fc = 1: 

C/fe(t, t') = 2T^T5tt' + A(0) (fc = 1) (B17) 
Cfc(i,t') = ^(e-'='l*-*'l-e-'='(*+*')) 



Let us now indicate the RG equation obeyed by the 
correlation function. It is obtained from considering the 



Using the equation for dtRk one can also write the closed 
equation for Ck{t,t') as: 
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,1 



1 



(B19) 



r 



and initial condition ljB18|l . Alternatively, one can work tions, including the determination of the FD violation ra- 
with Ck and Vk, which have a more complicated equation tio X{t, t') in the various regimes is left for forthcoming 



but simpler initial conditions: 

Vk{t,t') = 2rjT6u' (fc = l) 
Cfc(t,t') = §(e-'=^l*-*'l-e- 



publications. 



(B20) 



One easily checks that upon the assumption of time 
translational invariance as should hold in the equilibrium 
regime, the equation for Ck{t,t') = Ck{t — t') becomes, 
as expected, equivalent to the one for Rk{t — t') via the to lowest order in T coming from the cross term FA 
FDT relation. Further study of the nonequilibrium equa- reads: 



APPENDIX C: CORRECTIONS TO F-TERM BY 
PINNING DISORDER 

Let us give some details about the calculation of the 
graphs in Figs. l3l4l The correction to the effective action 



(5r = (i I {iUrt + i5Urt){iUrt' + iSUrf) '^{Urt - Urf + (5(u,.f - Urf)) 

^ Jr,ri,t.,t' 

Jti Jti Jti 

with Zr = iiirtdtUrt and the averages over 5u, 6u are restricted to 1-particle irreducible graphs. This splits into 
contributions corresponding to graphs (a,b,c) in Fig. ^ which evaluate respectively as (dropping all terms which do 
not correct F): 



= -iA"(0)i?2^=0 fi4F"[Zr]- f iUrt'iUrtidtUrtf) (C2) 
^ Jr Jtt' 

= i / A{Urt -Urt') < iSUrtiSUrt' F[Zr^ + j iUr^tidtJUr^ti] > 

^ Jr,ri,t.t' Jti 

= \R^q.uj=Q j dtdt'{A{Urt - Urt'))iUrtiUrt'F"[Zr] (C3) 

^p(c) — _2 / iu^^ < iSUrt'A{Urt — Urt' + SUrt — SUrf) 

2 Jr,ri,t,t' 

F[Zri + / 5iUr,tidt,Ur,ti + / ^ti (^Mri ti «Writi ] >= ^q,i^=0 / iUrtiUrt' dtUrtOf A' {Urt - Urt')F" [Zr] (C4) 
Jti Jti Jr,t,t' 

as well as the graph in Fig. ^ 

SY(3) ^ ^2 f iUrt < iSUrt'A{Urt - Urt' + SUrt - SUrt')F[Zri + ( SiUntidtiSUnti] > 

^ Jr,r\,t,t' Jti 

iUrtdt'A'{Urt-Urt')F'[Zr] j Rq^t-tiRq,ti-t' (C5) 



This yields the result H48|l in the text. 



APPENDIX D: MAPPING OF RANDOM lems, such as, after disorder averaging, the statistical me- 

FRICTION MODEL ONTO POLYMER AND 
RELATED PROBLEMS 



The random friction model (in its non-trivial T — 
limit) can be mapped formally onto various other prob- 
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chanics of a pure self interacting chain (e.g. a self avoid- 
ing walk problem) or, prior to averaging, to some random 
diffusion models, e.g. depolarisation of a spin diffusing 
in a random magnetic field. Concerning its behaviour 
one should distinguish between the genuine model (with 
a fixed distribution P{ri)) and the effective one which ap- 
pear as a coarse grained version of the pinning problem, 
in which P{ri) flows and becomes very broad. 

First setting P{r,t) ^ r]{r)urt one sees that (^5]) (with 
f{r,u) — and T = 0) is the Fokker- Planck equation 
dtP = VD{r){V + VV{r))P for the diffusion of a par- 
ticle with a random diffusion coefficient D{r) = l/rj{r) 
in a random potential V{r) = —lnrj{r), of equilibrium 
measure Peq{r) — e^^'^'"-' — rjlr). When ri{r) is uncorre- 
lated from site to site one does not expect any anomalous 
behaviour in any dimension, except if the distribution of 
rj has broad tails (e.g. algebraic would yield anomalous 
power law diffusion). In the effective model V{r) be- 
comes gaussian and grows with scale which corresponds 
to a particle localized in some regions of space. 

A complementary picture can be developed based on 
a mapping onto a self-interacting chain. The response 
function Rrt,r't' — dSurt / dh\h,=o of this model is obtained 
by solving: 



{v{r)dt - \7^)6urt = hS{r - r')5{t ~ t') 



(Dl) 



with initial condition 5urt=o = 0. This implies 5urt = 
for all t < t'. Thus the response is a function oi t — t' 
alone and its Laplace-Fourier transform s — iuj in any 
given random environment, can be written as: 



Rrr'{s) —< r\ 



I 



— -|- sr]{r) 



r' > 



du < r|e~"-"|r' > 



(D2) 



where H = — + sri{r), which has a positive spec- 
trum for s > ~s*. The value s* at which H develops 
an eigenstate of zero eigenvalue (e.g. s* = l/rjmax in 
the "classical" limit c 0) gives the large time decay of 

We can also write, in the Fourier domain, using the 
Feynman Kac formula: 

/'+OC px{u)—r 

Rrr'ittu)^ diie-(^+*'^^)" / Dx{v) 

Jo Jx(0)=r' 

eM-^ dv{-{-f + tu;Sv{x{vm (D3) 



with the "time" variables u and v. We have splitted 
ri{x) ~ r] + Sri{x) for convenience and added a small mass 
term /i for convenience. In this form the problem has the 
form of a spin decoherence problem, the integral being 
dominated by paths which average well over the random 
relaxation times, rather than paths with multiple returns 
to the same region which average poorly and then can- 
cel incoherently (details about the mapping and special 
distribution of noise can be found in^I). 



Averaging over disorder leads, for small disorder, to: 

Rr-r'.^ - / due-(^+'"'')«Z(r - r', u,g^ -uj^m) 
Jo ^ 

px{u)—r 

Z{r-r',u,g)^ / Dx{v) (D4) 

Jx{Q)=r' 

1 dx 

exp(— / dv-{ — )'^— / dvdv' g5[x{v) — xiv'))) 
Jo 4 dv 7o 

which is the partition function of a self avoiding walk 
in the Edwards representation. We have retained only 
the second moment 772 of Srjix), but of course the full 
interaction could be written using F[z\. The theory is 
described by a non-trivial fixed point in d < 4 (related 
to the n = 0{n) model with mass -\- iufff) and cou- 
pling g ~ f^A^^*^). This is compatible with the previous 
conclusions for the F theories using perturbation theory 
if we take t as a simple index, with no power counting 
dimension (it plays a role somewhat similar to the replica 
index a) . The correction to g by g^ comes from the con- 
traction of two rj2 vertices, of the form 5{Ld'^rj2) = ^^^2 
and is indeed logarithmically divergent in = 4. Note 
that the mass term is thus relevant at the fixed point 
9 = 9*- 

This analysis yields information at finite time. The 
Ginsburg criterion gives the critical regime as r > 
g-i/(4-d) or u > (^^2/(4-(i)^ which Z(r — r', M, g) takes 



the scaling form 



Z{r,u,g) = u 
Z{q = 0,u,g) 



-ud 



F[ru-'']Z{q^Q,u,g) 



(D5) 
(D6) 



with Sc{g) ~ eg for d > 2, Sc{g) = cg\n{l/g) for d ~ 2 
and Sc{g) = cg^^^ for d = 1. Thus for rf > 2 we expect 
that: 



1 



9''^ ~ — I 2 I 2 



in the non-critical regime, while we expect: 



Rq=0 



q=0,uj — 1 — 2 

' trjuj + crj2Uj 



(D7) 



(D8) 



i.e, Rq=o,t ~ p-'-e~'^^^*''^ in the critical regime, and for 
g > 0, the appropriate scaling function of rt^"^ . 
The critical regime correspond to: 



V2 



(D9) 



which for the genuine model gives a singularity only at 
finite (true) time (see however— for possibly more radical 
effect of non-gaussian disorder) . However in the limit of 
very broad disorder, as in the effective model, one has 
772 » rf and thus the critical singularity moves to small 

UJ. 
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Note that for d < 2 the behaviour is more radical as 
one expects, e.g. in d = 1: 



Rn 



1 



(DIG) 



Thus to conclude, in the absence of pinning disorder 
at T = the F term is preserved but generate higher 
order time derivative terms. The theory can be rescaled 
so as to possess a non-trivial finite w, finite disorder term, 
which presumably in d > 2 produces only preexponential 
algebraic corrections to the leading behaviours given by 
the most relevant rj term. For the effective model this 
critical behaviour should be observable even at small uj. 



APPENDIX E: FULL FLOW OF THE 
TRUNCATED EFFECTIVE ACTION 

In this appendix we discuss an approach to the cal- 
culation of the mean response function extending the 
approximate FRG scheme of Section III, but still ne- 
glecting the functional dependence of operators consid- 
ered in Section IV. In Section III the broad distribution 
of time scales was embodied in the so-called F term. 
We uncover here an interesting structure of additional 
operators in the spirit of the more complete set of mo- 
ments (i^Oi ■ ■ ' {t^")L discussed in the introduction. Re- 
call that in Sec. IIII Bl we showed that, although the F 
term did not renormalize itself, it did generate higher 
derivative terms. Such terms can contribute to the fre- 
quency dependence of the response function. We studied 
in the previous Appendix the response function of the 
pure random friction model, which does generate higher 
derivative terms, but neglects the scale dependence gen- 
erated by the pinning disorder. Here we consider these 
two effects in tandem, hence modifying the results for 
Rk{t), RkiLo) within the purely random friction model. 
While we have not been able to obtain simple expres- 
sions in this more complete (albeit still non-functional) 
approximation, one can go quite far in reducing the prob- 
lem to one of applied mathematics. 



1. Generalized random friction model 

To proceed further one needs to construct a more sys- 
tematic approach where all possible important terms in 
the dynamical action functional are included. We will 
generalize the F term (at zero temperature for simplic- 
ity) in the form: 



'S'kin — 



E E 

n— 1 pi ■■■'Pn>l 



F, 



(«) 



lUrti ■ ■ ■ lUrt„ 



(El) 



and we will often denote by m = pi -(- • • • -I- p„ the total 
number of time derivative in a given term of the sum. 



n=4 
m=7 



FIG. 12; Graphical representation of an F vertex Fp^ p^ 
with n = 4, pi = 1, p2 = 2, ps = 1, p4 = 3, m = J2iPi = 
Dots represent the number of time derivatives (i.e power of 
frequency factor), each leg has a different frequency. 



O 




FIG. 13: Compact notation for a generic F vertex. The 
open circle represents an F vertex with an arbitrary number 
of legs n not shown. On incoming (w) lines, an arbitrary 
number of time derivatives (powers of ui) are indicated by an 
open square. 



This form of the action neglects terms with products of 
time derivatives of Urt at the same time, as well as sta- 
tistically translationally invariant functional dependence, 
e.g. on Urii —Urt2- However, it does include considerably 
more physics, and as we will see, enough generality to ap- 
proach the problem of computing the averaged response 
functions. This kinetic part of the action corresponds to 
the following generalization of the random friction model: 



E 



+ F{urt,r) 



(E2) 



with [VpAn)-VpArn)]c = n!(-l)"+iF« . 

The response function is related to the lowest (n = 1) 
member of this hierarchy via 



OC 



(E3) 
(E4) 



and within the Wilson scheme the true physical response 
function R'i^^{iuj) is obtained via the same formula using 

the running |/=i,i(A/fc)- This is represented graphi- 
cally in Fig. [131 

One can carry perturbation theory using the general- 
ized F. The vertices are shown in Figs. I12I13I In this 
notation, there are a variety of important one loop di- 
agrams to be considered. These are shown in Fig. 1151 
Schematically, these contributions give rise to an RG 
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(a) 



FIG. 14: Graphs with only n = 1 vertices which enter in the 
response function (a) and its self energy (b). 
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(b) 



(c) 
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FIG. 15: One-loop diagrammatic contributions to the F 
terms. The diagrams in (a)-(e) represent contributions to 
F of order F, AF, F^, F^A, and F^A"^, respectively. 



equation for the Fm ^ of the form 



m' m—m' 



2 pC"') pC"-"') 



-(•E5) 



In (|E5p . we have neglected coefficients, powers of I, and 
fine distinctions such as the precise form of A which ap- 
pears in a given term. Repeated (primed) indices other 
than m and n are summed. The pure random friction 
model contains only those diagrammatic contributions 
with no pinning disorder A = 0. Thus only the second 
and third terms in (|E5|) are taken into account for in- 
stance by the mapping to a polymer problem in the pre- 
vious appendix. Note that there are no terms of O(FA^) 
contributing to the renormalization of F, which follows 
diagrammatically because one needs at least two F's to 
put boxes (time derivatives) on the incoming legs origi- 
nating from the two pinning disorder vertices. 

A beautiful simplicity arises due to the extremely 
broad distribution of timescales in the RG, and can be 
seen from the structure of ljE5|l . In particular, we note 
that the total number of powers of lo is conserved by 



all terms. Moreover, from the prior analysis, we expect 



7] 



(m) 



We thus conjecture that all F, 

p(n) _ „2m' 



in) 



scale 

in this way independently of n, i.e. Fm' ~ rji"^ ~™ (see 
below to see that this scaling is indeed self-consistent). 
Under this assumption, the terms involving more than 
two F's in l|E5|l can be seen to be strongly subdominant, 
which follows from the convexity of the 2m^ — m factor 
in the exponential. Thus the super-exponential (Gaus- 
sian) growth of the moments of the timescales, which is 
directly connected to the broad distribution of relaxation 
times, plays a key role in simplifying the structure of the 
RG. 

One can thus restrict the analysis to the linear part in 
F of the full one-loop RG equation. Note that this indeed 
combines the effects of pinning disorder and the "up- 
ward" feedback of the random friction model - the first 
two terms in (jESp . The linearized RG equation reads: 



(E6) 



diF^^l.p^ = ri{2n'-n)F^'^l.p^ 

n Pr — 1 

+ (n+l)ae-('^-2)'V y , , 

r—1 q—1 

where a = A'^^^A^. The general study of this equation is 
again highly difficult but we see that it does have special 
solutions where the F^^p^y depend only upon rii = X^iP*- 
From the above consideration about the asymptotic be- 
haviour it is rather natural to look for such solutions. 



Thus we let F, 
E6I) becomes: 



" l,m 



Then the hnear RG equation 



= r,(2n2-n)F(::] + (n+l)«e-(''-2)'(m-n)^^/;+i\ 



■(E7) 

This is much easier to solve and the attentive reader will 
easily find that this infinite hierarchy of differential flow 
equations is solved asymptotically by the ansatz 



^Pl-Pn ~ 



P_^(d-2+0)l 

a 



2m —m 



(E8) 



where m = "Yl^^iPi- The a^^ coefficients are given by 
ml^l 1 m!(2m + 2n-3)!! 



2m + 2r - 1 



n!(4TO- 3)!! 



The f^o™^ coefficients are not determined by the asymp- 
totic analysis. In principle, they should be matched at 

(n) 

some scale /* > to the form of the Fm coefficients de- 
termined by the early stages of renormalization, in which 
the linearized RG equation used to obtain them is not 
valid. One might imagine beginning with a model in 
which the bare relaxation time was distributed with cu- 
mulants f7o™^ ^'^'i naively ?7o™^ « We will use this 
prescription below purely in order to simplify notation. 
However, a more detailed analysis of the early stages of 
renormalization is in fact required to discern whether this 
is indeed correct. 
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2. Equation for the equilibrium response function 

We are now in a position to extract the response func- 
tion. Applying ljE4IE8IE9|) ior n — 1 gives the response 
as an infinite series in iuj: 



(2m2 - m)Uk 



+ [7- {d-2 + e)l-ln{l/iuj)]m 



,(E10) 



where Uk = C//=in(A/fc)- 

Since it is not easy to resum this series, and its con- 
vergence properties are unclear, we now reformulate the 
above calculation in a functional way, in hopes of sur- 
mounting the limitations of the expansion in terms of 
moments of relaxation times. Let us introduce the gen- 
erating function 

G{tuj,z)^ ^ M™-"^"F^„ (Ell) 

m > n > 1 

which conveniently captures both the F term (and hence 
distribution of relaxation times) and the response func- 
tion: 

F{z) = G(0,z), (E12) 
d 

+ Sfe(zcj) = — G{tuj,z)\ r.. (E13) 
oz 



Multiplying ljE7p by the appropriate powers of z and iuj 
and summing gives the flow equation 

diG = Ti{zd,G + 2z^dlG) (E14) 

This is somewhat simplified by defining the derivative 

Hiiu), z) — dzG{iLo, z), 

diH = Ti{H+5zd,M+2z^dlH)+ae-^'^-^'^'iLudi^ {iujd.H) 

(E15) 

A hopefully illuminating change of variables is to define 

u = \n{l/iuj), (E16) 
V = \n{z/iLu)-{e + d-2)l±\n{a/x)- (E17) 

Then K{u, v) — H{iuj, z), and obeys 

du,K ^ K+3dyK + 2d^,K+{du-dy) [e-'^d^K) , (E18) 

where diUi — xe^' as before. Note that the F term is 
recovered in the limit v ^ 00 and u — v = \nz + const, 
is fixed, in which the second term is negligible. In that 
limit, we recover the diffusion with drift equation H51|l . 
and the appropriate solution is K{u,v) = $(u — u) = 
F'(e"~"). More general solutions of l|E18() remain to be 
found. 



APPENDIX F: ONE LOOP HIERARCHY: METHOD OF CALCULATION 



In this Appendix we show how the systematic calculation of the one loop correction to the dynamical effective 
action can be organized, and sketch explicit calculation on the simplest examples. We focus on T = 0. 

The schematic form of the dynamical action S is given in the text in (|77|) as a sum of terms containing an increasing 
number of independent times (cumulants): 

S = iuiSi - -iuiiu2Si2 - -iuiiu2iu3Si23 - ■■ (Fl) 
2 D 

We use the same schematic notation where the indices 1, 2, 3.. are short hand notations for ^1,^2, t^.., space coordinate 
and all time and space integrations are implicit. From H77() Si is parameterized by an infinite set of kinetic coefficients 
77, D, .., S12 by a set of second cuniulant functions A, G, A, B, C, .., 5i23 by a set of third cumulants H, W, .. etc.. (from 

In a first stage we write the total one loop corrections to the action as the sum of tadpoles, two vertex loop, three 
vertex triangles etc.. with either iSi2 or 1S123 (and so on..) type vertices using the full response function R12, inverse of 
iuiSi, to contract the vertices (internal lines). Enumerating possible contractions and performing some combinatorics, 
yields upon grouping resulting terms by number of independent times: 

SSi = -< iu2Si2 > (F2) 

SS12 =< »lt3'5l23 > +[^'^34 < Suiusiu^ > -t- < 524^3 >< SisiUi >] (F3) 

1 3 

^^123 = [-^45 < Si23iU4iU5 > -f 3 < 514^5 >< 'S'235iW4 > +2 5i2ZU4«U5 > 5345] (F4) 

+ [2 < S34iu5 >< 5i5TO6 >< 526*^4 > +3 < Si2iuiiue >< S34iu5 > S^e] 

where additional (time) indices are integrated over Note that Si is corrected only by tadpoles, S12 by tadpoles and 
two vertex loops, and so on.. 



33 



In this formula notatfons such as, e.g. < Si2U3U4 > denote the sum of all possible contractions of the u fields with 
the u fields inside the bracket (at T = these are the only possible contractions). Since S12 — Si2[ui2, ui,U2,ui,U2, ■■] 
is an explicit function of U12 — ui — U2, and time derivatives of ui and U2 (and similarly for all other vertices) one 
may write the sum of all possible contractions as: 

< iU2Si2 >- (i?i25„,, + {dxRx2)du, + {dlRi2)du,) S12 + (F5) 

We recall that here causality R22 — restricts contractions only with ui (first term), iii (second term) etc.. 
The next stage is to make apparent rj, D etc.. and thus to define the expansion: 

+00 

Ri2 = -^{6i2 + Y,Aj,dl5,2) (F6) 

,2 +00 +00 

Ai = , A2 = -k-^D + k-'^rf (F8) 

The momentum structure of the one loop diagrams being trivial, within Wilson one can replace k = Ai everywhere. As 
explained in the text this expansion in power of frequency can be done consistently and corresponds diagrammatically 
to expansion in number of dots. 

One then evaluate the contractions shifting time integrations. Let us illustrate this on simple examples. 

The corrections to ij and D can be obtained from IIF2p using (|F5|) . One has: 



SSi = -(Ri2du,, + {diRi2)du, + {dlRi2)du,)Si2 (F9) 
= -[(<5i2 + A^di5i2 + A2dl5i2)du,,Si2 + {di5i2 + Aidl5i2)du,Si2 + dl5i2du^Si2] (FIO) 
= -,5i2[(l + Aid2 + A2dl)du,,Si2 + {d2 + A^dl)du,Si2 + dldu,Si2] (Fll) 



In the second line we have used the expansion (|F6|I and in the last line we have used that diRu = —d2Ri2 and 
integrated by parts over t2 ■ Then time derivatives on the vertex iSi2 can be evaluated and tranformed into derivatives 
w.r.t. fields as: 

82 = 82812 = {-U2dui2 + U29iiJ5i2 (F12) 
dl = a|5i2 = {-U2du,, + uldlJSi2 (F13) 

At all stages of the calculation we can drop all terms containing more than a fixed number (here 2) of time derivatives 
since they will contribute only to higher order terms in the effective action. Putting everything together we obtain: 

Sr]^G'{Q)-7]A"{0) (F14) 
6D = -A{0) + C"(0) - 277fc-2G"(0) - A"iO){k-^D - fc- V) (F15) 

which, in terms of rescaled quantities, yield the (|90I96|I in the text. 

Next we want to evaluate the corrections to second cumulant functions A, G, A, B, C encoded in SSu- 
We start by the simplest, the tadpole, which yields the feedback of third cumulants into second ones. 

5\tadpoleSl2 —< S123U3 >= (^13i9mi S'123 + diRi^du^Sus + -R235ti2 'S'l23 + 92-R239m2 S'l23 (F16) 

= 25i3[(l + Aid3)du, + dA,]Si23 (F17) 
= 2Si3[dui + Aiiisdusdui + U35«3^«i]'S'i23 (F18) 

with ^3 — u^du-i- This gives the 25*^(0,0, u) term in the equation (|68() for A and the H feeding term in the equation 
for G (we have not explicitly computed the feeding of W into A, S, C but it is easily obtained from the above). 

Next we need the 5^2 corrections to S12 which yield all non linear terms in H91I97I98I99|) for G,A,B,C. The 
corresponding correction to SS12 consists in the two terms in the square bracket in (|F3I) . The full calculation being 
tedious we only indicate here how one shuffles time integrals in the first term (denoted 6iSi2). Starting from: 

S1S12 = ^834 < S12U3U4 ^534[(i?14 " i?24)9„i2 + dlRl4du, + d2R24du2 + dfRlidu, + 9|i?249iiJ (F19) 

X[(i?i3 - i?23)9„i2 + (9ii?i3aui + 92i?239u2 + (9ii?139ui + i?239i.2]5i2 
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using again identities such that diRu = ^d^iRu, expanding the R and integrating by parts over and ^3 one can 
write 

S1S12 = ^534[(1 + Aidl + A2cf/){Sii - S24)dui2 +04+ ^iSi^) ((^M^ii + S24du2) + {SuOu, + (5249iJ] 

x[(l + Ai% + A2%''){di3 - S23)du,2 + {% + Aids^KSisdu, + ^2^ + %\hzdu^ + <5235fiJ]5i2 (F20) 

This is then in the form where, as above, aU time derivatives can be replaced by derivatives over fields acting either 
on S12 or 534 using identities such as (|F12p together with 9493534 = —u-iU/id^,^^Szi- The evaluation of the second 
term in ljF3p proceeds similarly and the sum of the two yields the (I91I97I98I99|) in the text. Note that causality must 
be enforced at each step of the calculation. 



APPENDIX G: INTEGRABLE 
"UNIRELAXATIONAL MODEL" 

In this section we introduce a set of integrable models 
in various dimensions which can be used as a check of the 
FRG equations derived in this paper. This models have 
a remarkable property that despite being random the dy- 
namics is extremely simple and the relaxation time scales 
are simply that of a single mode generalized oscillator. 

Let us consider first the toy model in zero dimension: 

ri{ut)dtut = f{ut) - rn^ut (Gl) 
with /(u) = ~V'{u), when: 

i^{u) = 77(1 - m- V'N) (G2) 
it can be rewritten: 

777 

dt{f{ut) - m'ut) = - — {f{ut) - m^ut) (G3) 

■n 

which can be integrated exactly, yielding a pure expo- 
nential relaxation with a single time scale: 

f{ut) - rn?ut = e"*^ (/(^o) - m^uo) (G4) 

Drawing F{u) = f{u) — rn^u as a function of u we see 
that all initial conditions starting in an interval between 
two adjacent minima and maxima and which contains a 
zero of F{u), will converge exponentially to this zero of 
Thus, if F{u) has several zeroes the convergence 
will be to minima and maxima of the potential energy 
H{u) = V{u) -\- ^m^u^ (depending on initial condition). 
This should not be a surprise since in that case 'q{u) 
changes sign so the dynamics is no more dissipative. One 
can then extend the system to non-zero temperature T > 
0, imposing FDT with a stationary measure e"^*^"^/^ by 
adding C,t to the right hand side of (|Gip with correlations: 

(CtCf) = 2r77(u) (G5) 

Note that this means imaginary noise at points where 
r](u) is negative. Thus the convergence to the maxima of 
V{u) is killed by interference effects. 

The exact response function associated to the u field 
can be obtained for this model at T = 0. Upon adding 



an infinitesimal perturbation ht on the r.h.s. of (|Gip the 
change 5u is such that: 

where Ut is given by (inH) and RfJ ^ ie"*"^. The 
disorder averaged response function is thus: 

i?,, = (l-/>^)-ii?(°) (G7) 

In the large time limit time translational invariancc is 
restored and one finds: 

Rw = R[°J (G8) 
This a consequence of the following property: 

(l-^^)-i = l (G9) 

where for each realization of the random fimction /(u), 
Ut is the solution of: 

fiut) = m^ut (GIO) 

and the average is taken w.r.t. any translationally invari- 
ant distribution for /(u). 

A similar model may be introduced in arbitrary dimen- 
sion. It is defined as: 

r]Frt = -Frt + Crt (Gil) 
i^rt = VVt + /(«rt,r) (G12) 

This yields the equation of motion: 

-rjV'^Urt - r]f'{Urt, r)Urt = V'^Urt + f{Urt, r) + Crt 

(G13) 

This is identical in form to the model discussed in the 
text with the exception that the damping coefficient is 
wave vector q dependent and vanishes as q^. Upon aver- 
aging over disorder one obtains an MSR action identical 
to the model studied in the text apart from the q^ mean 
damping with: 

G{u)=r]A'{u) (G14) 
A{u) = -fA"{u) (G15) 
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and no other higher order vertex for a gaussian dis- 
tributed / (more general expressions can be easily ob- 
tained for non-gaussian distributions). The bare response 
function of this model factors as: 



1 



+ 7]iuj) 



(G16) 



Similar arguments as above yield that this is also the 
exact response. 

The one loop Wilson FRG of this model is very similar 
to the one performed for the model in the text. Since the 
relaxation time is dimensionless in this model vertices 
such as G and A scale identically to A. Hence the appro- 
priate rescaled functions for these vertices are G ^ A;^*^G, 
A ^ ^r'^A. The one loop FRG equations are identical 
to the one given in the text for G and A apart from the 
(linear) rescaling part (not involving Q being identical 
to that for A. One can then check that the relation: 



G(u) = A'{u) 



(G17) 



A{u) = -lfK"{u) 



(G18) 



specific to this model, are indeed exactly preserved by 
the FRG, as announced in the text. Computing the cor- 
rection to the self energy yields: 



9,I](c.)-A,Af-2[A"(0)(i?(u;: 
+G'{Q){2iujR{uj) - iLuR{0)) + 



l-i?(0)) (G19) 
A{0)uj^R{lu)] 



One then checks that this exactly vanishes using A{0) = 
-^2A"(0), G"(0) = ^A"(0) and the above exact form for 
R{u;). 

This model can be further generalized to include sec- 
ond time derivative terms D ^ 0. Adding the term DF^t 
to the l.h.s. of (|nTT|) one obtains the model in the text 
with rj — > q^rj and D — > q^D in the bare inverse re- 
sponse function. Similar arguments yield invariance of 
the FRG function within the manifold IjlUUI ) given in 
the text (third cumulants have also been included). 
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